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Introduction 



The rich phenomenology of weak decays has always been a source of information 
about the structure of elementary particle interactions. A long time ago, (3- and \x- 
decay experiments revealed the nature of the effective flavour-changing interactions 
at low momentum transfer. Today, we are in a similar situation: Weak decays of 
hadrons containing heavy quarks are employed for tests of the standard model and 
measurements of its parameters. In particular, they offer the most direct way to de- 
termine the weak mixing angles and to test the unitarity of the Kobayashi-Maskawa 
matrix. On the other hand, hadronic weak decays also serve as a probe of that part 
of strong-interaction phenomenology which is least understood: the confinement of 
quarks and gluons into hadrons. In fact, it is this intricate interplay between weak 
and strong interactions that makes this field challenging and attractive to many 
theorists. 

Over the last decade or so, a lot of information has been collected about heavy 
quark decays from experiments on the T(4s) resonance, and more recently at e + e~ 
and hadron colliders. This has led to a rather detailed knowledge of the flavour 
sector of the standard model and many of the parameters associated with it. There 
have been several great discoveries in this field, such as B^-Bd mixingSH, b — > u 
transitions! - !, and rare decays induced by penguin operators!. Yet there is much 
more to come. Hopefully, the approval of the first 5-meson factory at SLAC has 
opened the way for a bright future for S-physics. At the same time, upgrades of 
the existing facilities at Cornell, Fermilab, and LEP will provide a wealth of data 
within the coming years. 

However, my main purpose here is to talk about theory, and fortunately there has 
been a lot of progress and enthusiasm in this field in recent years. This is related to 
the discovery ofheavy quark symmetr and the development of the heavy quark 
effective theory^ - ", which is a low-energy effective theory that describes the strong 
interactions of a heavy quark with light quarks and gluons. The excitement about 
these developments is caused by the fact that they allow (some) model-independent 
predictions in an area in which "progress" in theory often meant nothing more than 
the construction of a new model, which could be used to estimate some strong- 
interaction hadronic matrix elements. Therefore, I hope that you do not mind it 
if the part about the heavy quark effective theory constitutes the main body of 
these notes; this is where the main progress has been achieved from the theoretical 
point of view. Also, I have recently finished a long review article on heavy quark 
symmetryt3, which clearly simplifies the task of writing up lecture notes. 

In the first lecture, I will discuss the theory of heavy quark masses, their defi- 
nition in perturbation theory, and their determination from QCD sum rules. The 
second lecture provides an introduction to the standard model description of quark 
mixing, different parametrizations of the Kobayashi-Maskawa matrix, the status of 
the determination of the mixing angles, and the physics of the unitarity triangle. 
The third lecture, which covers more than half of these notes, is devoted to a re- 
view of the ideas on heavy quark symmetry and the formalism of the heavy quark 
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effective theory. In particular, I will discuss the theory of the model-independent 
determination of | V c b\ from exclusive semileptonic B — > D*£ v decays. At the end of 
each lecture you will find some suggestions for little exercises, which are typically 
related to the derivation of important equations given in the notes. I have tried to 
select problems that are fun to solve and contain some interesting pieces of physics. 
You are invited to see if I am right. 
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1. Heavy Quark Masses 



Because of confinement at large distance scales, quarks and gluons do not appear 
among the physical states of QCD. There is thus no natural, physical definition of 
quark masses. Rather, several definitions are possible and have been proposed, and 
it is often a matter of convenience which one to use. Most of these definitions are 
tied to the framework of perturbation theory. In this lecture, I will discuss some of 
the most common mass definitions and their interrelation. Special emphasis is put 
on the discussion of the running quark mass in the context of the renormalization 
group. I will then discuss how values for the bottom and charm quark masses are 
obtained from QCD sum rules. Very briefly, I will touch upon the subtle issue of 
an infrared renomalon in the pole mass, which implies an intrinsic uncertainty in 
these determinations. 

1.1. Quark Mass Definitions 

For heavy quarks, the nonrelativistic bound-state picture suggests the notion of 
a pole mass rnQ le defined, order by order in perturbation theory, as the position of 
the singularity in the renormalized quark propagator. One can show that the pole 
mass is gauge-invariant, infrared finite, and renormalization-scheme independent^. 
It is thus a meaningful "physical" parameter as long as the heavy quark is not 
exactly on-shell. For instance, the pole mass appears in the formula for the energy 
levels in quarkonium systems. 

An alternative gauge-invariant definition is provided by the running quark mass 
in some subtraction scheme. One usually works with the modified minimal subtrac- 
tion (MS) scheme^ and denotes the running mass by mg(yu). Running quark masses 
are used when there is a large momentum scale fi ^> tuq in the problem, so as to 
absorb large logarithms, which would otherwise render perturbation theory invalid. 
Of course, one can use other subtraction schemes such as minimal subtraction^! 
(MS). Sometimes, it may even be convenient to use a gauge-dependent definition 
of a heavy quark mass such as the so-called Euclidean mass mg ucl , which is defined 
by a subtraction at the Euclidean pointi! p 2 = —rrvq. What is important is that 
these perturbative definitions are related to each other in a calculable way, order 
by order in perturbation theory. 

One may also think of obvious nonperturbative definitions of a heavy quark mass. 
For instance, one could define tuq to be one half of the mass of the ground state in 
the corresponding (QQ) quarkonium system, or as the mass of the lightest hadron 
containing the heavy quark, or as that mass minus some fixed binding energy, etc. 
The problem is, of course, that these definitions are ad hoc, and relating them to 
each other or to the perturbative definitions given above would require to solve 
QCD, a task that is presently beyond our calculational skills. There is, however, a 
bridge between the two classes of definitions, which is provided by the heavy quark 
effective theory, which is an effective field theory appropriate to describe the soft 
interactions of an almost on-shell heavy quark with light quarks and gluons. There, 
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one can define a parameter A, which corresponds to the effective mass of the light 
degrees of freedom in a hadron Hq containing a single heavy quark Q, in terms of a 
gauge-invariant^infrared finite, and renormalization-scheme independent hadronic 
matrix element". One can then define a heavy quark mass m*q as^ll uiq = mn Q — A. 
In perturbation theory, and up to corrections of order I/tuq, the mass tuq coincides 
with the pole mass. However, the above definition is more general, as it does not 
rely on a perturbative expansion. 



1.2. Quark Masses in Perturbation Theory 

Let me now discuss the concept of heavy quark masses to one-loop order in 
perturbation theory. The bare quark mass mg are appearing in the QCD Lagrangian 
is related to the renormalized mass mq n by a counter term, 

m h c? e = - Sm Q , (1) 

where mg re and 5itiq are divergent quantities. The counter term is chosen such 
that the renormalized mass is finite. Its value depends, however, on the subtraction 
prescription. In the MS scheme, the renormalized mass will depend on some sub- 
traction scale /i, i.e. m™ = mg(/i). To calculate Smq in perturbation theory, one 
has to evaluate the heavy quark self-energy E(p) shown in Fig. [I]. The relation is 

5niQ = divergent part of E(j/ = mq) + scheme-dependent finite terms. (2) 

To calculate theself-energy at one-loop order, it is convenient to use dimensional 
regularizationE^fi, i.e. to work in D = 4 — 2e space-time dimensions. Then the 
renormalized coupling constant is related to the bare one byEj 

a!? are = Z a a a (ji) fi 2€ = a s fi 2e + 0(a 2 s ) . (3) 

From a straightforward calculation, one obtains at one-loop order 

£y/= m o) — m Q tt~ ( 2 ) r ( e ' 



3vr \A7Cfi 2 J w 1 - 2e ' 
^^f- + ln 4 + ^+0(e), (4) 



7r \e m , 



where 1/e = 1/e — 7 + ln47r. Notice that this result is still /i-independent, since 
the /1 dependence of a s /e cancels the /^-dependent logarithm. A scale dependence 
appears only when one subtracts the ultraviolet divergence. In the MS scheme, the 
mass counter term is defined to subtract the 1/e-pole in the self-energy, i.e. 

m Q (» = m Q (l + ^ . (5) 
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The pole mass, on the other hand, is defined so as to absorb the entire self-energy 
on-shell. Hence 

5m p Q ole = E(#= m Q ) , 



m^ = m Q l + ^fJ + ln^ + Df. (6) 



7r \e mg 3 
Comparing the two results, one obtains 

r n »oie f 1 4a a (m Q ) ) 

m Q\ m Q) - m Q Y g~ j> 

m Q (» =m Q (M Q ) 1 1- — In ^-l. (7) 

The first equation relates different perturbative definitions of the heavy quark mass, 
namely the pole mass to the running mass in the MS scheme evaluated at the scale 
jj, = mQ. I have written this relation as a perturbative expansion in powers of 
OL s (jnn). This expansion is also known to two-loop order from a calculation by Gray 
et al£3 Similar relations exist in other subtraction schemes, e.g. 

MS, MSN pole/, OsKf)/4 \<l 

m Q( m Q)= m Q I 1 ^^V3 _7+ Jj' 

EucU Euch pole d a s{ m Q ) n i 

™Q ( m Q ) = m Q I 1 ^ 21n2 

where the gauge-dependent Euclidean mass is defined in the Landau gauge 




em = -+mmm+- = — ^ L — + 



Figure 1: The quark self-energy to one-loop order in QCD. 

The second relation in (0) determines the running of the heavy quark mass. In 
QCD, quarks become "lighter" at high-energy scales. The simple one-loop calcu- 
lation presented above is, however, not sufficient to scale tUq{^jl) up to very large 
scales. For instance, using (|7|) to extrapolate the mass of the bottom quark up to a 
typical grand unification scale Mqut ~ 10 16 GeV, one would obtain a meaningless 
result * 

m^rrib) it ml n 

Of course, the problem of large logarithms in perturbation theory is not specific to 
this case, but is encountered frequently. Large logarithms can be controlled in a 
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systematic way using the beautiful machinery of the renormalization group. It is 
worth going through the solution in great detail. 



1.3. Renormalization- Group Improvement 

For the purpose of this section, I will switch to a slightly more general notation 
and rewrite (17f) as 



m(fi) = m(/i ) (l - — In ^ + 0(a 2 s ) } . (10) 

I 7T fi Q J 

If n ^> /io, the logarithm can be so large that a^ln/i 2 //! 2 , becomes of order unity, 
and ordinary perturbation theory breaks down. In fact, one can show that at any 
given order in perturbation theory there will be large logarithms of the type 

(j?y^ ( In n 2 /fi 2 y 



where I have used the fact that the running coupling constant scales like a s (fi) ~ 
1/ ln(/i 2 /Aq CD ). It is necessary to resum these "leading logarithms" to all orders in 
perturbation theory. This is achieved by solving the renormalization-group equation 
(RGE) for the running quark mass, 

li m(//) = 7(a s ) m(//) . (12) 

d/i 

The anomalous dimension 7 has a perturbative expansion in the renormalized cou- 
pling constant: 



7(«.) =70^ + 7! +.-.. (13) 

The coefficients 7, are known to three-loop orderllEi For our purposes, however, 
it is sufficient to note that0 

404 40 
7o = -8, 71 = — 3~ + ~9~^' 

where rif denotes the number of quark flavours with mass below //. Throughout 
these notes I will evaluate the QCD coefficients for N c = 3 colours, and only display 
the dependence on the number of flavours explicitly. The value of 70 follows from 
the one-loop result (|T0|). 

The next step is to rewrite the RGE in the form of a partial differential equation, 
making explicit the scale dependence of the renormalized coupling constant. This 
gives 

/ d d \ 

V~djl + ^ a '* da (jjS ~ 7( ~ a ' V m ^ = ° ' ^ 15 ^ 
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The (3- function 



(5(a s 



f ■ 



djjL 



-2a, 



a-.. 



47T/ 



+ 



(16) 



describes the running of the coupling constant. The one- and two-loop coefficients 
arel§ _ 

A) = ll-^n/, ^ = 102-^/1/. (17) 
The exact solution of the RGE can be written in the form 

m{fj) = U(n,fj,o)m(no) , (18) 
with the evolution operatoi@ilIl 



U (fjL, fj, ) = exp J da 



7(q) 
(3(a) 



(19) 



The trick is to obtain a perturbative expansion in the exponent of this expression 
by inserting the expansions for the anomalous dimension and /3-function. After a 
simple calculation, one finds 



a a {lJL) J 



1 + 



a B (/x ) - a s (yu) 7i A) - #i7o 



47T 



2/5 2 



(20) 



In this result, the running coupling constant has two-loop accuracy. The corre- 
sponding expression is 



oi s (jj) 



Air 



A)hi(/i 2 /A 2 



QCDJ 



1 - 



ft lnlnQu 2 /A| CD ) 
/3 2 ln(^/A| CD ) 



(21) 



where Aqcd is a scheme-dependent scale parameter! The factor containing the ratio 
of the running coupling constants in (|20| ) sums the leading logarithms (|TT| ) to all 
orders in perturbation theory. Keeping just this factor corresponds to the so-called 
leading logarithmic approximation (LLA), which is often used to attack problems 
containing widely separated scales. Note that to calculate this factor, one only needs 
to compute the one-loop coefficient of the anomalous dimension. This is a rather 
trivial task, as it suffices to calculate the 1/e-pole in the quark self-energy. Once 
the leading scaling behaviour has been factored out, perturbation theory becomes 
well-behaved, i.e. the next-to-leading corrections in the parenthesis in (|2U|) obey 
a perturbative expansion free of large logarithms. This is why the approach is 



^The value of a s (no) at some reference scale n$ can be used to eliminate the dependence on Aqcd- 
Nowadays, it is convenient to choose = mz, since a s (mz) is known with high accuracy. 



7 



called "renormalization-group improved perturbation theory" . The terms of order 
a s , which I have shown explicitly, correspond to the next-to-leading logarithmic 
approximation (NLLA). Including them, one sums logarithms of the type 



a s [a s In (22) 

to all orders. To achieve such an accuracy, it is necessary to calculate the two-loop 
coefficients of the anomalous dimension and /5-function. 

Combining the above results and setting fiQ = mg, one obtains the running 
quark mass to next-to-leading order in renormalization-group improved perturba- 
tion theory. The result is 



4/A) 



».,W^W h + S ^Sst^M +0{ M, (23) 



7T 



where 



S = — + — ( 24 ) 



5 34 107 

It is to be understood that the number of flavours changes as fi crosses a quark 
threshold^. For instance, when one uses ( ]2"B| ) to scale the running bottom quark 
mass up to large scales, nf changes from 5 to 6 when /i crosses the top quark mass. 



At the same time, the QCD scale parameter Aqcd in expression ( |21| ) for the running 
coupling constant changes, so that a s (fi) is a continuous function. From the point 
of view of convergence of perturbation theory, eq. (^3j) can be used to evaluate the 
running mass at an arbitrarily large scale fi. However, I have derived this result 
assuming that there are only QCD interactions. But at high-energy scales, other 
gauge and Yukawa interactions become important as well. They are most significant 
for the running of the top quark mass, as I will discuss in detail in the next section. 
In the case of the bottom quark, the result that renormalization effects lower the 
running mass at high-energy scales remains true. In fact, in many extensions of the 
standard model it is possible to obtain a unification of the bottom quark and the 
tau lepton masses at a typical grand unification scale: 

m b (M GUT ) = m r (M GUT ) for M GUT ~ 10 16 GeV . (25) 

For more details on this, I refer to the lectures by L. Hall in this volume. 



1-4- Running Top Quark Mass 

The interplay of gauge and Yukawa interactions leads to very interesting effects 
on the evolution of the top quark massElEll. In my discussion below, I will focus 
on the standard model. The analysis for extensions of the standard model such as 
supersymmetry proceeds in a similar way&~lli 
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Let me write the running top quark mass in terms of the Yukawa coupling At(//): 



m t (M)=M M )^, (26) 

where v denotes the vacuum expectation value of the Higgs field, normalized so that 
v(mw) — 246 GeV. It is convenient to define a running coupling constant «*(//) as 

««(A0 = ^ • ( 2 7) 

The most important contributions to the evolution of the running mass m t (fi) come 
from the QCD interactions as well as from the large top Yukawa coupling. The 
relevant one-loop diagrams are depicted in Fig. ||. They lead to the RGE 

d , . . ct q (u) 3 OiAu) 

"d^ km ' M = - 8 i^ + 2^ + - (28) 

where the ellipses represent other contributions, which I will neglect. Combining 
this with the RGE for v(n), 

^^ = -3^ + ..., (29) 

one obtains 

/x — InaJ/i) = -16 — h9— h . . . . (30) 

d/i 47T 47T 

To get an idea of the structure of this equation, suppose first that the QCD coupling 
constant is not running, i.e. a s = const. Then the RGE has an infrared- stable fixed 
point ata f = y« s , meaning that irrespective of the value of at(/i) at large scales, the 
Yukawa coupling is attracted into the fixed point as \i decreases. This is illustrated 
in Fig. |3|. The resulting fixed-point value of the top quark mass, m t = ^2na s v, 
is entirely determined by the low-energy group structure of the theory, which is 
responsible for the one-loop coefficients in (|3Q|) . 



H 



+ — — 1 >— — + 



Figure 2: Dominant one- loop contributions to the self-energy 
of the top quark in the standard model. 
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ln(\)JM) 

Figure 3: The fixed-point structure of the RGE (^) in the 
case a s = const., with 4a s /7r = 0.1. 



For running a s (fi), one has to find the simultaneous solution of ( |30| ) and the 
evolution equation for the strong coupling constant, 

,A. lnaM = .U^l. (31) 

Depending on the initial conditions, one finds that there are still "quasi-fixed point" 
solutions, where the running of at(^) is connected to the running of a s (/i), but inde- 
pendent of the value of the Yukawa coupling at large scales. In leading logarithmic 
approximation, the exact solution is not too hard to obtain. It reads 



8/7 



a t (M) \a s {M) 



1 9 g t (M) 
2 a s (M) 



MM), 



1/7 



(32) 



where M 3> fi denotes some large mass scale, at which the initial conditions are 
imposed. In order to illuminate the structure of this equation, let me distinguish 
three cases: 



(i) a t (M) < a s (M): 

In this case, the top Yukawa coupling is weak, and 



a t (fi) ~ a t (M) 



a s (jd) 



8/7 



a a {M)J 

which is nothing but the standard QCD evolution [cf. 



can be approximated by 

(33) 
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(n)a t (M) = la s (M): 

This case corresponds to the quasi-fixed point obtained by Pendleton and RosseS, 
where (1321) simplifies to 

2 

«*0) = g as{n) ■ (34) 
However, to obtain this solution requires a fine-tuning at the large scale M. 



0.3 




1 -^^^ 

0.02 0.04 0.06 0.08 0.1 0.12 

a s ([L) 



Figure 4: Evolution of the coupling constants for different 
initial values of ctt(M). I assume that M = 10 16 GeV. 



(iii) a t (M) > a s (M): 

This is the most interesting case of a large top Yukawa coupling. Then the approx- 
imate solution of (1321) reads 



s a a (n) _ 



1/7' 



(35) 



This is the quasi-fixed point discovered by HillEl. At low scales fi <ti M, a t (fi) 
follows a s (fi) and becomes independent of the initial value at(M). The dynamics 
that generates the top Yukawa coupling at some large scale M cannot be tested. An 
example of this mechanism is provided by the top condensation model of Bardeen, 
Hill, and Lindner^. Using (|35|), one obtains for the fixed-point value of the top 
quark mass 



m t {m t ) 



2v(m t ) 



yV a s (m t ) 



1 



'qs(M) 
a s (m t ) 



l/7-i -1/2 



210 GeV . 



(36) 



where I have assumed M ~ 10 16 GeV. Again, the result is almost entirely determined 
by the group structure of the low-energy theory. The dependence on M is very weak. 
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In Fig. |], I show the evolution of the coupling constants according to ( |32| ) for 
a set of initial values of a t (M). The quasi-fixed point behaviour is clearly visible 
once a t (M) is large enough. A more careful analysis including all standard model 
contributions in the RGE givesc3 m t ~ 225 GeV (for M ~ 10 16 GeV). Somewhat 
smaller values m t < 200 GeV are obtained in the minimal supersymmetric standard 
modelE3~E2l, depending however upon the value of tan (5. 



1.5. Determination of and m c from QCD Sum Rules 

In this section, I will discuss the extraction of the masses of the bottom and 
charm quarks from QCD sum rules, and in particular from the analysis of quarko- 
nium spectra. The main idea of sumrules was developed in the pioneering papers of 
Shifman, Vainshtein, and ZakharovEHl. Their idea was to use quark-hadron duality 
to obtain a prediction of hadron properties from calculations in the quark-gluon 
theory of QCD. Consider, as an example, the vacuum polarization induced by the 
electromagnetic current of a bottom quark: 

% J d 4 x e** ( | T {r(x),f(0)} | ) = (q^ - g^q 2 ) U(Q 2 ) , (37) 

where j M = 67^6, and Q 2 = —q 2 . The Lorentz structure of the correlator follows 
from current conservation. The function H(Q 2 ) satisfies a once-subtracted disper- 
sion relation 

-d^ nw2) ^/ ds (jrW Imn(s) ' (38) 

where by the optical theorem 

13s 

is related to a measurable cross section. Here, e^ = —1/3 is the electric charge of 
the 6-quark, and "6 b" is a short-hand notation for "hadrons containing b 6" . In the 
sum rule analysis one considers the moments of the correlation function, which are 
given by 

1 



M, 



n(Q2) lQ^o = ^/ dsS ^ llmII(s) - (40) 



n! V dQ 2 , 

In principle, these moments can be extracted from experiment. 

As long as Q 2 is not close to the resonance region, the function H(Q 2 ) can be 
calculated in perturbative QCD. Since Q 2 = is far away from the physical threshold 
for bottomonium production, this assumption is justified in the present case. Hence, 
perturbative QCD should provide a good approximation to a calculation of the 
moments M n . The leading and next-to- leading perturbative contributions to the 
correlator are shown in Fig. [5[ They give rise to the spectral density^ 



1 v(3 -v 2 ) ^. 2 . r 4a 



Imn pert ( S ) = — v - ' &(v 2 ) 1 + 



4tt 2 I 3tt 



vr 2 3 + v /7T 2 3 



2v 4 V 2 4 



(41) 
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where v — y 1 — Am 2 /sis the relative velocity of the heavy quarks. In this expression 
one uses the Euclidean mass for in order to minimize the effect of radiative 
corrections [see (El)]. For s ^> 4mf, the perturbative spectral density leads to the 
well-known cross section 

, i , T . s»W Ana 2 9 / a s \ , 

a s (e+e- ^66) pcrt - 6 — iV c e^l + -i). (42) 

However, the above form of the spectral density includes threshold effects, which 
become important when s comes closer to the threshold region. 



O <D <2> <S> 

Figure 5: Perturbative contributions to the correlator U(Q 2 ). 
The current operators are represented by circles. 

In QCD sum rules, one usually adds to the perturbative contributions nonper- 
turbative corrections, which are power suppressed and appear at higher orders in 
the operator product expansion (OPE) of correlation functions such as U(Q 2 ). In 
the case at hand, the leading nonperturbative corrections are proportional to the 
gluon condensate (a s G 2 ) ~ 0.04 GeV 4 , which in sum rules is treated as a phe- 
nomenological paramete ri. In the case of the bottomonium system, one is in the 
fortunate situation that the nonperturbative contributions to the spectral function 
are very small, as they are suppressed by the ratio (a s G 2 ) jm\ ~ 10~ 4 . Their effect 
on the mass of the bottom quark is below 1% and can safely be neglected. 

The idea is now to construct two equivalent representations for the moments 
M n , a "theoretical" one based on expression ( fj"T| ) for the spectral function, and a 
"phenomeno logical" one based on a measurement of the e + e~ — > bb cross section: 

1 °° 

M* h ~ - / d SS - n - 1 Imn pcrt ( S ), 

7T J 

Am 2 b 
1 °° 

AC P = - [ ds s^- 1 Imn cxp ( S ) . (43) 

7T J 

The moments M* h are very sensitive functions of m^. By comparing them to a 
phenomenological expression that uses detailed experimental information in the b b 
vector channel, one can extract an accurate value for the bottom quark mass. Exper- 
imentally, six resonances have been identified in this channel. Their masses Mr and 
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electronic widths r R (e + e ) are known rather precisely^. To evaluate the integrals 
over the experimental cross section, one then makes the following approximation: 

Imn exp ( S ) ~ -4-2 E M R r R (e + e-) 5(s - M R ) + ImII pert ( S ) &(s - s ) , (44) 

where the sum is over the narrow resonances T(ls) to T(6s), and the continuum 
above some threshold value so is approximated by perturbation theory. This can 
be justified by quark-hadron duality. Inserting this ansatz into (|43|) and equating 
the two representations for the moments, one obtains the sum rules 

So 

- J ds S — 1 Imn pert ( S ) ~ — — 2 £ M^ 1 r R (e + e") . (45) 



TT < R 

The goal is to find an optimal set of parameters (mj, s ) such that M* h ~ M^ xp for 
many values of n. In practice, the perturbative calculation breaks down for large 
values of n, since the Coulomb corrections proportional to a s /v in (|4l]) then become 
too large. Typically, one is limited to values n < 10. As mentioned above, the non- 
perturbative corrections are very small (below 1%) in the case of the bottomonium 
system. They become sizeable, however, in the charmonium system, where a simi- 
lar analysis can be performed. Therefore, mj can be extracted with larger accuracy 
than m c . 

Let me then quote the numerical results obtained by various authors following 
the strategy outlined above. Values_fpr the Euclidean mass of the bottom quark 
have been obtained by Shifman et al£^l hjif ucl = 4.23 ±0.05 GeV), Guberina et al.E3 
( m Euci = 4.i9±o.06 GeV), and Reindersll (mf ud = 4.17±0.02 GeV). The differences 
in these values are mainly due to changes in the experimental input numbers. The 
result of Reinders is the most up-to-date one. Values for the Euclidean mass of 
the charm quark have been obtained many years ago by Novikov et al£3 (m^ ucl ~ 
1.25 GeV) and Shifman et alB (m^ ucl = 1.26 ± 0.02 GeV). They have not changed 
since then. Converting these results into pole masses using (|8|), one findU 



m P° le = 4.55 ±0.05 GeV 



mP° le = 1.45 ± 0.05 GeV . (46) 

The increase in the error is due to an additional uncertainty in the value of a s (mg). 

Let me finally mention some alternative approaches to extract pole masses from 
QCD sum rules. Voloshin has proposed to resum the large Coulomb corrections 
proportional to (a s /v) n to all orders in perturbation theory, using a nonrelativistic 
approachil. This allows one to go to higher moments, thereby suppressing the 
resonance contributions to the sum rule. The disadvantage is that some relativistic 
corrections are not taken into account. The result is a rather large value of the pole 
mass of the bottom quarkEZl, 

ml° lc = 4.79 ± 0.03 GeV . (47) 
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Another alternative is to obtain heavy quark masses from the study of sum rules 
for heavy-light hound states such as the B- and 5*-mesons. From such an analysis, 
Narison obtainsE2 

m P° le = 4.56 ±0.05 GeV. (48) 
A more recent analysis using the heavy quark expansion to order 1 / tuq givesSSS 

m b = 4.71 ±0.07 GeV, 

m c = 1.37 ±0.12 GeV. (49) 

It should be noticed, however, that these sum rules are more sensitive to nonper- 
turbative corrections than are the sum rules for the quarkonium systems. 

1.6. Infrared Renomalons 

As emphasized at the beginning of this lecture, the fact that at low energies 
quarks are always confined into hadrons prohibits a physical, on-shell definition 
of quark masses. Although for heavy quarks the notion of a pole mass is widely 
used, such a concept becomes meaningless beyond perturbation theory. No precise 
definition of the pole mass can be given once nonperturbative effects are taken into 
account. 

It is interesting that indications for an intrinsic ambiguity in the pole mass can 
already be found within the context of perturbation theory, when one studies the 
asymptotic behaviour of the perturbative series for the quark seif-eiiergyEJl53. It 
can be shown that the existence of so-called infrared renomalonsE3~E3 generates a 
factorial divergence in the expansion coefficients. Roughly speaking, the reason 
is that self-energy corrections to the gluon propagator in the diagram depicted in 
Fig. |l] effectively introduce the running coupling constant a s (Vk 2 ~), where k 2 is the 
square of the virtual gluon momentum in Euclidean space. Since a s (^/k 2 ~) increases 
as k 2 decreases, small loop momenta become more important. One may reorganize 
the perturbative expansion as an expansion in the small coupling constant a s {mQ) 
using 

^-".WEI^W^ • ( 5 °) 



n=0 

„2 



When one then performs the loop integral over k , the extra logarithms give rise to 
a growth of the expansion coefficients proportional to n\. In such a situation, it is in 
principle not possible to improve the accuracy of perturbation theory by including 
more and more terms in the perturbative series. Starting from some order n, the 
size of the corrections increases. The best that can be achieved is to truncate the 
series at an optimal value of n. This introduces an intrinsic error, which depends 
exponentially on l/a s (mg). In the case of the pole mass, one can show that the 
irreducible uncertainty in the value of Wq o1c is of orderl3i3 

a <° = 4 mQ exp ( - a«) * ■ (5i> 
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It is thus of the order of the confinement scale Aqcd- The fact that ^.wlq Q / rriq e 
vanishes in the tuq — > oo limit justifies the concept of a pole mass for heavy quarks, 
at least in an approximate sense that is limited to the context of perturbation 
theory. 

Although one has to be careful when using (^Tj) to obtain an estimate of the 
intrinsic uncertainty in the pole mass, I will insert Aqcd ~ 150 MeV to find Amg° le ~ 
50 MeV. This uncertainty should be kept in mind when considering numerical results 
for heavy quark masses obtained, for instance, using QCD sum rules. 

1.7. Exercises 



Derive the one-loop expression (Q) for the on-shell quark self-energy in QCD. 
The necessary loop integrals in D space-time dimensions can be found in any 
reasonable textbook on quantum field theory. 



Show that ([19]) is the exact solution of the RGE (fTpD, and derive the next-to- 
leading logarithmic approximation ([2Up for the evolution operator. 

By calculating the 1 / e-poles of the diagrams shown in Fig. |^, derive the one- 
loop RGE (ESI) for the running top quark mass. 



Solve the RGE (|30|) for the case a s = const., and show that there is an 
infrared-stable fixed point at a t = ^§a s . 



Derive eqs. (HHg5|). 
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2. Quark Mixing in the Standard Model 

In this lecture, I give an introduction to flavour-changing decays and quark mix- 
ing in the standard model. I will discuss quark mixing in the cases of two and three 
generations, introduce some useful parametrizations of the Kobayashi-Maskawa ma- 
trix, and briefly review the status of the direct experimental determination of the 
entries in this matrix from tree-level processes. I will then turn to the geometrical 
interpretation of the mixing matrix provided by the unitarity triangle, and finish 
with some remarks on rare decays. Since you probably have heard many lectures 
on these subjects, my presentation will be rather short. Let me refer those of you 
who want more details to two excellent review articles that I like very much: A Top 
Quark Story by Buras and Harlander£3, and CP Violation by Nirc^l. 



2.1. Cabibbo-Kobayashi-Maskawa Matrix 



Let me briefly remind you of some facts about the flavour sector of the standard 
model. Below mass scales of order m w ~ 80 GeV, the standard model gauge group 
SUc(3) x SUl(2) x Uy(l) is spontaneously broken to SUc(3) x U cm (l), since the 
scalar Higgs doublet acquires a vacuum expectation value 



(0) 



v [0 



246 GeV . 



(52) 



This gives masses to the W- and Z-bosons, as well as to the quarks and leptons. The 
quark masses arise from the quark Yukawa couplings to the Higgs doublet, which 
in the unbroken theory are assumed to be of the most general form that is invariant 
under local gauge transformations. The Yukawa interactions are written in terms 
of the weak eigenstates q' of the quark fields, which have definite transformation 
properties under SUl(2) x Uy(l). After the symmetry breaking, one redefines the 
quark fields so as to obtain the mass terms in the canonical form. This has an 
interesting effect on the form of the flavour-changing charged-current interactions. 
In the weak basis, these interactions have the form 



C - y 

*-int — 1= 



9_ 
V2 



fd'r 



W\ + h.c. 



(53) 



In terms of the mass eigenstates q, however, this becomes 



Ant = -— 7= c l, t L ) Y V } 



KM 




Wt + h.c. 

fj, 



(54) 



The Kobayashi-Maskawa mixing matrix 

^ V u d V us Vub ^ 

V K m ^ Vcd V cs V cb 

\V td V u Vtb) 



(55) 
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is a unitary matrix in flavour space. In the general case of n quark generations, 
Vkm would be an n x n matrix. I will now discuss the structure of this matrix for 
the cases of two and three generations. 



2.1.1. Mixing Matrix for Two Generations 

In this case, V is a 2 x 2 unitary matrix and can be parametrized by one angle 
and three phases: 



V 



cos 9 C e ia sin 9 C e l/3 \ 
-sm9 c e^ cose c e l(/3+ ^ ) 

e ia W cos#c sin# c \ ( 1 
e' 7 I 1 - sin 9 C cos 9 C / [ 



(56) 



The three phases are not observable, however, as they can be absorbed into a 
redefinition of the phases of the quark fields ul, cl, and sl relative to d^. After this 
redefinition, the matrix takes the standard form due to CabibbdZi 

V C =( 008 " c S[nd ff c ) . (57) 
V — sin C cos 9c J 

The Cabibbo angle 9c can be extracted from K 7re + z/ e decay. Experimentally, 
one findst3 sin^c — 0.22. 



2.1.2. Mixing Matrix for Three Generations 

A 3 x 3 unitary matrix can be parametrized by three Euler angles and six phases, 
five of which can be removed by adjusting the relative phases of the left-handed 
quark fields. Hence, three angles 9^ and one observable phase 5 remain in the 
quark mixing matrix, as was first pointed out by Kobayashi and MaskawaB. For 
completeness, I note that in the general case of n generations, it is easy to show that 
there are |n(n— 1) angles and — l)(n— 2) observable phasesO. Whereas therefore 
the original Cabibbo matrix was real and had only one parameter, the Kobayashi- 
Maskawa matrix of the standard model is complex and can be parametrized by four 
parameters. 

The imaginary part of the mixing matrix is necessary to describe CP violation 
in the standard model. In general, CP is violated in flavour-changing decays if there 
is no degeneracy of any two quark masses, and if the quantity J CP ^ 0, where 

J CP = | ImiVijVuVJVZj) | ; i±k, j^l (58) 

is invariant under phase redefinitions of the quark fields. One can show@ that all 
CP-violating amplitudes in the standard model are proportional to Jcp. 
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I will now discuss two of the most convenient parametrizations of the mixing ma- 
tig:. The "standard parametrization" recommended by the Particle Data Grouped 

/ C12 C13 S12C13 Si3e~ lS 



Vv 



KM 



-S\2 C23 — c 12 s 23 s 13 e%5 C12 C23 — S12 S23 S13 e l(5 S23 C13 
V s 12 ^23 — c 12 C23 S13 e* — C12 S23 — S12 C23 S13 e* C23C13 




Here, one uses the short-hand notation cy = cos 0y and = sinfly. 
tages of this parametrization are the following: 

• I V u b\ = S13 is given by a single angle, which experimentally turns out to be 
very small. 

• Because of this, several other entries are given by single angles to an accuracy 
of better than four digits. They are: V u d ~ c 12 , V us ~ si 2 , V c b — S23, and 
V* 

• The CP-violating phase 5 appears together with the small parameter S13, 
making explicit the fact that CP violation in the standard model is a small 
effect. Indeed, one finds 

JCP = I «13 S 2 3 S12 S S cl 3 C 23 C12 I • (60) 

For many purposes and applications, it is more convenient to use an approximate 
parametrization of the Kobayashi-Maskawa matrix, which makes explicit the strong 
hierarchy that is observed experimentally. Setting c 13 = 1 (experimentally, it is 
known that C13 > 0.99998) and neglecting s i3 compared to terms of order unity, one 
finds 

C12 S12 si 3 e" 

Vkm - I S12 C23 C12 c 23 s 2 3 j • (61) 

S\2 $23 — C12 C23 S13 e lS — C12 S23 C23 

Now denote si 2 = A ~ 0.22. Experiments indicate that S23 = 0(A 2 ) and S13 = 0(A 3 ). 
Hence, it is natural to define S23 = AA 2 and Si3e~* 5 = A\ 3 (p — irj), with A, p 
and 77 of order unity. An expansion in powers of A then leads to the Wolfenstein 
parametrizationEl 



V K m - 



I 1-f A AX 3 (p~i7])\ 

-A 1-f AX 2 

\ A A 3 (l -p-irj) -AX 2 1 / 



0(A 4 ) . (62) 



It nicely exhibits the hierarchy of the mixing matrix: The entries in the diagonal 
are close to unity, V us and Vd are of order 20%, V c b and V ts are of order 4%, and 
V u b and Vtd are of order 1% and thus the smallest entries in the matrix. Some care 
has to be taken when one wants to calculate the quantity Jcp in the Wolfenstein 
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parametrization, since the result is of order A 6 and thus beyond the accuracy of the 
approximation. However, taking i — u, j — d, k — t, and / = b in ([58]), one obtains 
the correct answer 

J CP ~ A 2 r]\ 6 ~ 1.1 x 10- 4 A 2 r], (63) 

which shows that Jqp is generically of order 10~ 4 for A ~ 0.22. As a consequence, 
CP violation in the standard model is a small effect. 

In principle, the elements in the first two rows of the mixing matrix are accessible 
in so-called direct (tree-level) processes, i.e. in weak decays of hadrons containing 
the corresponding quarks. In practice, as I will discuss below, the entries | V u d\ and 
| V^ s | are known to an accuracy of better than 1% from such decays, whereas | Vcd\, 
\V CS \, and | Vcb\ are known to about 10-20%. The element | V u b\ has an uncertainty of 
about a factor of 2. The same is true for | V t d\, which is obtained from a rare process, 
namely E>d-Bd mixing. There is no direct information on \V ts \ and \V t b\- The large 
uncertainty in | V u b\ and | V t d\ translates into an uncertainty of the Wolfenstein 
parameters p and i]. A more precise determination of these parameters will be the 
challenge to experiments and theory over the next decade. 

2.2. Experimental Information on Vkm from Tree-Level Processes 

In this section, I will briefly review what is known about entries of the Kobayashi- 
Maskawa matrix from tree-level processes. Rare decays, which are induced by loop- 
diagrams involving heavy particles, will be discussed later. 

2.2.1. Determination of | V u d\ 

From a comparison of the very accurate measurements of super-allowed_nuclear 
/3-decay with /i-decay, including a detailed analysis of radiative corrections^^, one 
obtains^ 

| V ud \ = 0.9744 ±0.0010. (64) 

2.2.2. Determination of | V us \ 

An analysis of K e % decays, i.e. K + — > -K°e + v e and K\ — > 7r~e + z/ e , including 
isospin- and SU(3)-breaking effects§, yields \V US \ = 0.2196 ± 0.0023. Alternatively, 
from an analysis of hyperon decays, including SU(3)-breakingcorrectionsE3, one 
obtains | V us \ = 0.222 ± 0.003. This leads to the combined valueH 

| K s | = 0.2205 ±0.0018. (65) 

Let me add a note on K e3 decays in this context. The fact that the theoretical 
uncertainty in the description of these exclusive hadronic decays is only about 1% 
is at first sight surprising. The reasons are the following: The approximate SU(3) 
flavour symmetry of QCD implies that, at zero momentum transfer, the hadronic 
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form factor f+^ n that parametrizes K —>■ it transitions induced by a vector current 
equals unity, up to symmetry-breaking corrections. The Ademollo-Gatto theorem£9 
states that these corrections are of second order in the symmetry-breaking param- 
eter (m s — m q ), where q = u or d. Chiral perturbation theory can be employed to 
calculate the leading symmetry-breaking corrections (chiral logarithms) in a model- 
independent way". In my third lecture, I will discuss that there are analogues to 
all these ingredients in the case of semileptonic 5-decays. 

2.2.3. Determination of | V c d\ 

From the combination of data on single charm production in deep-inelastic 
neutrino-nucleon scattering with the semileptonic branching ratios of charmed me- 
sons, one deduces^ 

\V cd \ = 0.204 ±0.017. (66) 



2.2.4. Determination of | V^ s | 

Here the usefulness of deep-inelastic scattering data is limited, since the extrac- 
tion of \V CS \ would depend upon an assumption about the strange quark content 
of the nucleon. It is better to proceed in analogy to the determination of |V^ S |, 
i.e. to use the semileptonic D e3 decays D — > Ke + u e . The experimental data are 
compared to various model calculations of the decay width^irtl This procedure 
leads, however, to a significant amount of model dependence. The result isEa 

| = 1.00 ±0.20. (67) 

You may wonder why the theoretical description of the decay width is so un- 
certain. The reason is that, in general, our capabilities to calculate hadronic form 
factors in QCD are rather limited. An understanding of the nonperturbative con- 
fining forces would be necessary to make the connection between quark and hadron 
properties, which is a prerequisite for any calculation of hadronic matrix elements. 
In the case of K e % decays, a symmetry of QCD helps us eliminate most of the 
hadronic uncertainties. For D e3 decays, there is no such flavour symmetry. Hence, 
one has to rely on phenomenological models. Let me briefly mention some of them: 

• Isgur, Scora, Grinstein. and Wise (ISGW) have proposed a nonrelativistic con- 
stituent quark modeEl They solve the Schrodinger equation for a Coulomb 
plus linear potential to obtain meson wave functions. Weak decay form factors 
are calculated at maximum momentum transfer, where both mesons have a 
common rest frame, by computing overlap integrals of the wave functions of 
the initial and final meson. The nonrelativistic approximation breaks down if 

5max' 

• Bauer, Stech, and Wirbel (BSW) have constructed a relativistic mo 
which uses light-cone wave functions to calculate weak decay form factors 
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at zero momentum transfer. In order to extrapolate to q 2 > 0, they assume 
nearest-pole dominance. 

• Korner and Schuleiil (KS) use a variation of the BSW approach, in which the 
q 2 dependence of theform factors is adjusted according to asymptotic QCD 
power-counting rulesE3. 

• There are avariety of other models, for instance by SuzukiS, Altomari and 
WolfensteintZI, and Grinstein et al.E3, which are more or less closely related to 
one of the above. 

Disadvantages of these models are that their relation to the underlying theory of 
QCD is not obvious, that it is hard to obtain an estimate of their intrinsic uncer- 
tainty or range of applicability, that they rely on ad hoc assumptions (for instance, 
concerning the q 2 dependence of form factors), and that it is difficult to improve 
them in a systematic way. 

Another stream of research focuses on analytical or numerical approaches that 
bear a closer relation to field theory. QCD sum rulesll provide a fully relativistic 
approach based on QCD. They rely, however, on the assumption of quark-hadron 
duality and need some phenomenological input. Nonperturbative effects are mod- 
elled in a simple way by introducing few universal numbers, the so-called vacuum 
condensates. There have been extensive applications of sum rules to the study of 
meson weak decay form factors^ - ". These analyses can be done for all physical 
values of q 2 . However, in spite of the many successes of QCD sum rules it must be 
said that the only known approach to low-energy QCD that truly starts from first 
principles is lattice gauge theory. For more details on this, I refer to the lectures by 
P. Lepage in this volume. Let me just mention that, as far as weak decay form fac- 
tors are concerned, these computations are extremely complex and CPU-consuming. 
They are not yet competitive with (less rigorous) analytical approaches. 

2.2.5. Determination of | V u b\ 

Since b — > u transitions are strongly suppressed in nature, their discovery in 
inclusive B — > X u i v decays was one of the breakthroughs in I?-physics in recent 
yearsErt]. In order to account for CP violation in the standard model, it is necessary 
that Vub be different from zero; otherwise Jqp = 0, and the observed CP violation in 
the kaon systemE3 could not be explained by the Kobayashi-Maskawa mechanism. 

The first direct observation of a b — > u transition was provided by two fully 
reconstructed events, one B° — > Tr"tu~u decay and one B + — > ufi + v decay, re- 
ported by the ARGUS collaboration^. However, the number of exclusive charmless 
5-decays that have been observed so far is too small to obtain a reasonable deter- 
mination of | V u b\. A high luminosity S-factory would improve this situation. Given 
enough statistics, the idea is to obtain a model-independent measurement of | V u b\ 
by using heavy quark symmetry to relate the decay form factors in B — > X u I v 
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and D — > Xd i v transitions^, at the same value of the recoil energy of the light 
hadron X q . In the limit of infinite heavy quark masses (mj, m c — ► oo), the ratio of 
these form factors is fixed in a model-independent way. In the real world, however, 
there are nonperturbative power corrections to this limit. As a consequence, even 
in an ideal measurement the ratio | V u b/V c d\ can only be determined up to hadronic 
corrections of order 

'Aqcd Aqcd" 



1 + c 



+ 



(68) 



with a coefficient c of order unity. Model calculations^^ show that these correc- 
tions are of order 15% for X = n or p. I believe that it is fair to say that even 
at a high-luminosity 5-factory, the prospects for getting a determination of | V u b\ 
from exclusive decays, which is more reliable than that from inclusive decays, are 
rather limited. Nevertheless, it is certainly worth while to pursue this strategy as 
an alternative. 




Figure 6: Sketch of the lepton spectrum in inclusive semilep- 
tonic .B-decays. I assume | V u b/V c b\ — 0.08. The b — > u signal 
has been multiplied by a factor of 5 to be visible on this plot. 



The present determinations of | V u t\ 
momentum spectrum in inclusive B —. 



are based on measurements of the lepton 
X q ti> decays, where X q is any hadronic 
state containing a g-quark, with q = c or u. The expected signal is shown in Fig. ||. 
Over most of the kinematic region, the spectrum is dominated by b —>■ c transitions, 
which are strongly enhanced with respect to b — > u decays. The only place to 
observe charmless transitions is in a small window above the kinematic limit for 
B — > D £ v, but below the kinematic limit for B — » it Iv decays: 
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(69) 
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i.e. 2.31 GeV < p e < 2.64 GeV. Indeed, the ARGUS and CLEO collaborations 
have reported signals in this regionta - 0. An extraction of | V u t\ from these data is 
difficult, however, as the shape of the spectrum close to the kinematic endpoint is 
dominated by nonperturbative effects. This can be seen as follows. The conventional 
description of inclusive decays of hadrons containing a heavy quark starts from the 
free-quark decay model, in which the hadronic decay B — > X q i v is modelled by the 
quark decay b — > qtv. One usually calculates the decay distributions to ordeiipfcin 
perturbation theory by including the effects of real and virtual gluon emissionE2il _ ll22l. 
Computing then the lepton spectrum, one finds that the kinematic limit for b — > u 
decays is given by p™ ax = m^/2 ~ 2.35 GeV, where I neglect the w-quark mass 
and take ~ 4.7 GeV for the sake of argument. The point is that the region in 
which the b — > u signal is observed experimentally is almost not populated in the 
free-quark decay model. Nonperturbative bound-state effects, such as the motion of 
the 6-quark inside the 5-meson, are responsible for the population of the spectrum 
beyond the parton model endpoint. The way in which such effects are incorporated 
into phenomenological approaches^S is to a large extent model-dependent. Only 
very recently has there been progress towards a model-independent description of 
the distributions near the kinematic endpoint in the context of It has 

been shown that the leading nonperturbative effects can be parametrized in terms 
of a universal structure function, which describes the momentum distribution of the 
6-quark inside the I?-meson. The hope is that it will be possible to get an accurate 
prediction for this function using various theoretical approaches. I believe that, 
within a year or two, it should be possible to reduce the theoretical uncertainty in 
the analysis of the inclusive decay spectrum by a factor of 2 or so. 

Alternatively, one can describe the lepton spectrum close to the endpoint as the 
sum of contributions from semileptonic -B-decays into exclusive final states. These 
exclusive decays are then described using one of the various bound-state modelsEirE-S 
discussed above. Of course, there is again a strong model dependence associated 
with this approach. 

Because of these theoretical uncertainties, the current value of | V u b\ has a model 
dependence of almost a factor of 2. In 1991, the ARGUS collaboration obtained! 



I Vub/Vcb\ = 0.11 ± 0.012, using the approach of Altarelli et alH!-3, which is based on 
the parton model. Using instead bound-state models to describe the spectrum by a 
sum over several exclusive decay channels, values for | V u b/V c b\ between 0.10 and 0.22 
are obtained. However, the most recent data reported by the CLEO collaboration 
lead to significantly smaller values!: | V u b/V c b\ = 0.076 ± 0.008 using the approach 



of Altarelli et aLllrJ, and 0.05 < | V u b/V c b\ < 0.11 using bound-state models. Thus, I 



think a reasonable value to quote is 

0.08 ±0.03. (70) 



Vub 

V cb 
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2.2.6. Determination of I V r 



cb\ 



As for | V^b |, it can be extracted from both inclusive and exclusive .B-decays. 
I start with a discussion of inclusive decays. The idea is to compare the total 
semileptonic branching ratio to the parton model formula 



Br(B ^X q lu) = r B [ Vc | V cb f f + Vu \ V ub f] , (71) 



G%ml f ,„ „ Jm 2 c 

where r] c and 1] U contain the short-distance QCD corrections, and 

f(x) = 1 - 8x + 8x 3 - x 4 - 12x 2 lnx (72) 
is a phase-space correction due to the mass of the charm quark. The QCD correction 



factor for b — *> c decays ~ 0.87. The contribution from b — > u transitions 

is very small and, at the present level of accuracy, can be neglected. 

Let me note that it has recently become clear how to compute in a systematic 
way the nonperturbative bound-state corrections to the parton model resul t ( jTlf) . 
For the total inclusive decay rates, these corrections turn out to be of ordertilFEj 
(Aqcd / "W*) 2 '■ Numerically, they are very small and play only a minor role in the 
analysis^. 

The disadvantage of using inclusive decays to extract | V c b\ is that there are 
substantial theoretical uncertainties due to the fact that the total rate scales as 
the fifth power of the bottom quark mass. For instance, an uncertainty of only 
±150 MeV in the value of m& translates into an uncertainty of ±8% in the value of 



\Vcb\. For this reason, the last time the particle data groupie used inclusive decays 
to obtain a value for | V c b\ was in 1988. A compilation of more recent analyses of 
inclusive decay spectra has been given by StonelH3. He finds 

I K^l f-^^y /2 = 0.040 ± 0.005 , (73) 
VI. 5 ps/ 

where I have normalized the result to a value of the 5-meson lifetime that is close 

as 



to the world average t b q = (1.48 ± 0.10) ps reported last year by Danilov 
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well as to the most recent value tbo = (1.52 ± 0.13) ps obtained from an average of 



LEP resultsliM The error quoted in (|73|) does not take into account the theoretical 
uncertainty in the value of m^. 

Ultimately, a more precise determination must come from the analysis of ex- 
clusive decay modes. With the development of the heavy quark expansion, it has 
become clear that certain symmetries of QCD in the heavy quark limit can help to 
remove much of the model dependence from the theoretical description of exclusive 
semileptonic -B-decays. This will be explained in detail in my third lecture. In 
particular, for the decay mode B — > D*£ v the theoretical uncertainties in the mo- 
mentum distribution of the D*-meson close to the zero recoil limit can be controlled 
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at the level of a few per centliM The mos t, recent analysis of this spectrum by the 



CLEO collaboration leads to the valued 



V, 



tb 



1/2 



<hl[ ^—j = 0.039 ±0.006, (74) 

which I will use in the analysis below. In contrast to (|73|), here the uncertainty 
is dominated by the experimental systematic errors and, to a lesser extent, the 
statistical ones. The theoretical uncertainty is well below 10%. There is thus room 
for a significant improvement of the accuracy within the next few years, as more 
data become available. 

Finally, let me mention that by combining the above result with the value of 
| V us | given in fl6"5]) , one obtains 



A 



for the Wolfenstein parameter A. 



2.3. Unitarity Constraints 



Vr 



cb 



V 2 
* us 



0.80 ±0.12 



(75) 



The fact that the quark mixing matrix must be unitary can be employed to 
put limits on the elements that are not easily accessible to experiments. One has 
to assume, however, that there are only three quark generations. Unitarity then 
implies that 

1-IK/-IK/- (76) 



V I 2 
V t j\ 



Inserting here the known values of | V u j | and | V c j\ , one obtains 



0.003 < | V td \ 
0.032 <\V ts \ 
0.9991 < I Vtb\ 



< 0.016, 

< 0.048 , 

< 0.9994. 



(77) 



Note, in particular, that | Vtb\ is constrained to an accuracy of better than 10 3 . 
Unitarity also imposes a rather tight constraint on | V cs \. Using that 



1 - IK 



cd | 



v cb \ 2 



(7* 



one finds 



0.976 < | V c 

This should be compared to the value (|6 



< 0.981 . (79) 
extracted from a direct measurement. 



2.4- The Unitarity Triangle 



A simple but beautiful way to visualize the implications of unitarity is provided 
by the so-called unitarity triangleEH, which uses the fact that the unitarity equation 



VijV* k 
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can be represented as the equation of a closed triangleJn, the complex plane. There 
are six such triangles, all of which have the same area 
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1 r 

2 CP 



Under phase reparametrizations of the quark fields, the triangles change their ori- 
entation in the complex plane, but their shape remains unaffected. 

Most useful from the phenomenological point of view is the triangle relation 



v ud v: b + v cd v; h + v td v f 



cdVcb 



'tdV t b 







(82) 



since it contains the most poorly known entries in the Kobayashi-Maskawa matrix. 
It has been widely discussed in the literaturdHZI _ t2j. In the standard parametriza- 
tion, V cd V* b is real, and the unitarity triangle has the form shown in Fig. 0a. It is 
useful to rescale the triangle by dividing all sides by V cd V* b . The rescaled triangle is 
shown in Fig. |7|b. It has the coordinates (0,0), (1,0), and (p,r)), where p and r\ ap- 
pear in the Wolfenstein parametrization (|62|). Unitarity amounts to the statement 
that the triangle is closed, and CP is violated when the area of the triangle does 
not vanish, i.e. when all the angles are not zero. 
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Figure 7: (a) The unitarity triangle; (b) its rescaled form in 
the p-rj plane. The angle 7 coincides with the phase S of the 
standard parametrization (|59|) . 

To determine the shape of the triangle, one can aim for measurements of the two 
sides R b and R t , and of the angles a, /3, and 7. So far, experimental information is 
available only on the sides of the triangle. The current value of | V u b\ in (170) implies 



0.24 < R h < 0.53. 



13) 



As mentioned above, the uncertainty in this number is mainly a theoretical one, so 
any significant improvement must originate from new developments in the theory 
of inclusive -B-decays. 
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Let me now discuss the determination of | V t d\ and of the side R t of the unitarity 
triangle. The process of interest here is the mixing between B^- and 5^-mesons, 
which in the standard model is a rare process mediated through box diagrams with 
a virtual top quar k exchange. The ARGUS and CLEO collaborations have reported 
measurementsl 135 H 136 l of the mixing parameter Xd defined as Xd = Am b c1 Ts d , where 
ArriB denotes the mass difference between the mass eigenstates in the Bd-Bd system. 
The weighted average of their results is Xd = 0.68 ± 0.08. Combining this with the 
-B^-meson lifetime^, r Bi = (1.48 ±0.10) ps, one obtains Am B = (0.46 ±0.06) ps" 1 . 
Recently, direct measurements of Am# have been reported by some of the LEP 
experiments. The average value presented at the Moriond Conference isEi^ Am# = 
(0.519tg;c)5g) ps -1 . I combine these results to obtain 

Am B = (0.49 ± 0.04) ps" 1 . (84) 

The theoretical expression for Am# in the standard model is 

G 2 m 2 

Am * = - J hF ^QCD m B {B B f B ) S{m t /m w ) \ V td V t * b \ 2 , (85) 

67T Z 



where t?qcd — 0.55 contains the next-to-leading order QCD corrections^, and 

, |l39| , |l40 
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S(m t /m w ) is a function of the top quark massl 139 H 140 l. For 100 GeV < m t < 300 GeV, 
it can be approximated bylill S{m t /mw) — 0.784 (rrit/mw) 1 ' 52 - The product (Bb/b) 
parametrizes the hadronic matrix element of a local four-quark operator between 
5-meson states. Recently, there has been a lot of activity and improvement in 
theoretical calculations of the _B-meson decay constant fs, using both lattice gauge 
theory and QCD sum rule calculations. For a summary and discussion of the 
extensive literature on this subject, I refer to the review articles by Buras and 
Harlanderii and by myselill. Here I shall use the range 

B B /2 f B = (200 ± 40) MeV , (86) 
which covers most theoretical predictions. Solving (|5]) for | V t d\, one obtains 

Using the numbers given above, as well as m t = (170 ± 30) GeV, I find 

\V td \ = 0.0087 ±0.0023. (88) 

The corresponding range of values for R t is 

0.75 < Rt < 1.44. (89) 

The error in this value will hopefully be reduced soon, with a measurement of the 
top quark mass at Fermilab. Further improvements could come from a measurement 
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of fs in B + — > t + v t decays, or from the observation of B S ~B S mixing. In the latter 
case, one could use the relation 

x jl + SU(3)-breaking corrections j (90) 

to obtain a rather clean measurement of | V^|. However, in the standard model one 
expects large values x s ~ 10-30, which will be very hard to observe in near-future 
experiments. 
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Figure 8: Experimental constraints on the unitarity trian- 
gle in the p-rj plane. The region between the solid (dashed) 
circles is allowed by the measurement of R^ (Rt) discussed 
above. The dotted curves show the constraint following from 
the measurement of the e parameter in the kaon system. The 
shaded region shows the allowed range for the tip of the uni- 
tarity triangle. The base of the triangle has the coordinates 
(0,0) and (1,0). 

In Fig. |8|, I show the constraints imposed by the above results for Rj, and R t in 
the p-r\ plane. Another constraint can be obtained from the study of CP violation 
in ^-decays. 

2.5. CP Violation in K -Decays 

The experimental result on the parameter e that measures CP violation in K-K 
mixing implies that the unitarity triangle lies in the upper half p lan e pr ovided the 
so-called Bk parameter is positive. Most theoretical predictions^ TL^J fall in the 
range Bk = 0.65 ±0.20, supporting this assertion. The arising constraint in the p-T] 
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plane has the form of a hyperbola approximately given byt 



7/ 



P )A 



2 ( m t 



1.52 



+ Pc 



A 2 B 



K 



0.50 



(91) 
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where P c ~ 0.66 corresponds to the contributions from box diagrams containing 
charm quarks, and A = 0.80 ± 0.12 according to fl75|). I have included this bound 
in Fig. |8|. 

In principle, the measurement of the ratio Ke(e'/e) in the kaon systemE 
could provide a determination of r\ independent of p. In practice, however, the 
theoretical calculations^ - llH of this ratio are affected by large uncertainties, so 
that there currently is no useful bound to be derived. Further information on the 
parameters p and r\ could be extracted from rare if-decays, such as Kl — > 7r°z/z/, 
K + — > 7r + z/z/, and Kl — > p + p~. The corresponding branching ratios vary between 
10 -11 and 10 -10 . The experimental detection of such decays will be very hard. 



2. 6. Potential Improvements in the Determination of p and r\ 

The main conclusion to be drawn from Fig. |^ is that, given the present theoretical 
and experimental uncertainties in the analysis of charmless 5-decays, B—B mixing, 
and CP violation in the kaon system, there is still a large region in parameter 
space that is allowed for the Wolfenstein parameters p and rj. This has important 
implications. For instance, the allowed region for the angle (3 of the unitarity triangle 
(see Fig. [?D is 6.9° < /3 < 31.8°, corresponding to 

0.24 < sin 2(3 < 0.90. (92) 

Below I will discuss that the CP asymmetry in the decay B — > ip Kg, which is one of 
the favoured modes to search for CP violation at a future 5-factory, is proportional 
to sin 2/9. Obviously, the prospects for discovering CP violation with such a machine 
crucially depend on whether sin 2/3 is closer to the upper or lower bound in (p2|). 
A more reliable determination of the shape of the unitarity triangle is thus of the 
utmost importance. 

Let me briefly discuss what I think are the potential improvements in this de- 
termination in the near future. It is not unrealistic to assume the following: 

• The top quark will be discovered at Fermilab. 

• The uncertainty in the value of | V c b\ will be reduced to ±0.04, mainly through 
a better control of the experimental systematic errors. 

• Both theoretical and experimental improvements will reduce the uncertainty 
in the determination of | V u b\ by a factor of 2. 

• Improved theoretical calculations, in particular using lattice gauge theory, 

1/2 

will allow the determination of the nonperturbative parameters Bg fg and 
Bk with an accuracy of 10%. 
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Figure 9: Constraints on the unitarity triangle obtained as- 
suming a more precise determination of the various input pa- 
rameters, as explained in the text. 

For the purpose of illustration, I shall assume that m t = 170 GeV, | V cb \ = 0.039 ± 
0.004, | V ub /V cb \ = 0.080 ± 0.015, B^ 2 f B = (200 ± 20) MeV, and B K = 0.65 ± 0.07. 
This leads to R b = 0.38±0.08 and Rt = 1.09±0.20. Similarly, the constraint derived 
from the measurement of e becomes more restrictive. This fictitious scenario is 
illustrated in Fig. [| Obviously, the parameters of the triangle (the value of r], 
in particular) would be much more constrained than they are at present. Notice 
also that once a sufficient accuracy is obtained, there is a potential for tests of the 
standard model. If the three bands in Fig. [| did not overlap, this would be an 
indication of new physics. 

2.7. Rare B -Decays 

Let me briefly touch upon the interesting subject of rare -B-decays, i.e. de- 
cays induced by lo op d iagrams. There are nice review articles on this subject by 
Bertolinili^B and Alili^H. Rare -B-decays are dominated by short-distance physics, 
namely the top quark contribution in penguin and box diagrams. The decay rates 
are usually rather sensitive to the mass of the top quarkMll. Long-distance effects 
are much less important than in rare i^-decays. Hence, rare -B-decays are "clean" 
from the theoretical point of view. They are ideal to determine the parameters of 
the unitarity triangle. In addition, these decays can provide important tests of the 
standard model, since they are often sensitive probes of new physics. 

So far, the only rare .B-decay that has been observed isQ B — > K*j. In general, 
decays of the type B — > X s 7 are mediated by penguin diagrams such as the ones 
shown in Fig. |10|. These decays receive very large, but calculable, short-distance 




31 



QCD correct ionsll^3 _ E2j, which strongly enhance the decay rate. In addition, it is 
important to take into account the effect of gluonbremsstrahlung, which affects the 
total decay rate as well as the photon spectrum!^. All such effects are reasonably 
well under control for the inclusive B — > X s ^ decay rate. For m t = 170 ± 30 GeV, 
the prediction isEa 

Br(B -> X s 7) = (2.8-4.2) x 10~ 4 . (93) 

Of more interest from the experimental point of view are exclusive rare decay modes. 
Since B — > K^y is forbidden by angular momentum conservation, the lightest final 
state that can appear is K*j. To obtain from fl9~3|) a prediction for Br(B — > K*~f) 
requires an estimate of the ratio 

T(B -> A-7) 

R "--r(B^x,j)- (94) 

This is where hadronic uncertainties enter the analysis. Estimates of Rk* presented 
in the literature vary between 4% and 40%. I believe that a reasonable value ii^ni^a 
R K * ~ 10-15%. This leads to 

Br(B -> K* 7 ) ~ (3-6) x 10~ 5 . (95) 

This is in agreement with the result reported by the CLEO collaboration!: 

Bi(B -> K*j) = (4.5 ± 1.9 ± 0.9) x 10" 5 . (96) 




b 
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Figure 10: Penguin diagrams contributing to the rare decay 

Further rare -B-decays, for which there are "clean" theoretical predictions (which 
depend, however, on the top quark mass), include B — > X s v v with a branching ratio 
of 10~ 5 -10~ 4 , B — > X s £l with a branching ratio of about 10~ 5 , and B s — > ££ with 
a branching ratio of about 10 -7 . 

2.8. Measurements of the Angles of the Unitarity Triangle 

Let me finish this short introduction into the phenomenology of rare decays 
with a very beautiful application, namely the direct measurement of CP violation 
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in decays of neutral 5-mesons into CP eigenstates /. For these processes, there are 
often no (or only very small) hadronic uncertainties. The CP asymmetries arise 
from an interference of mixing and decay. When such a decay is dominated by a 
single weak decay amplitude, the time-integrated asymmetrv obeys a very simple 
relation to one of the angles of the unitarity triangleSl^i^: 



jdt \r(B°(t) - /) - r(B°(t) - /) 

1 

oo r _ 

jdt r(B°(t)-/) + r(B°(t)-/) 



± sin(2<y9) 



X 



1 + x 2 



(97) 



where x = Arris t~b- Examples of such decays are: 



B d 




(-sin 2(3) , 


B d 


^4>K S 
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(-sin 2/3) , 
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-> 7T + 7T~ 


(sin 2a) , 


B s 
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(— sin 27) , 


B s 


^ct>K s 


(sin 2(3) . 



(9? 



The decay mode B d ^ ip K$ is particularly "clean" and is often referred to as the 
"gold-plated" mode of a future £> -factory. 



2.9. Exercises 



• Show that the generalization of the Kobayashi-Maskawa matrix to n quark 
generations can be parametrized by \n{n — 1) angles and |(n — l)(n — 2) 
observable phases. 

• Show that the quantity Jqp defined in (|58| ) is invariant under phase redefini- 
tions of the quark fields. Verify that in the standard parametrization Jcp is 
given by flSOD . 

• Derive that the maximum lepton momentum in the decay B — > X £ v is given 
by \m B {l- m 2 x /m B ) [cf. (§9|)]. 

• Convince yourself that all six unitarity triangles have the same area (|8T]) . 

• Show that the decay B — > K'y is forbidden by angular momentum conserva- 
tion. 
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3. Heavy Quark Symmetry 



In this lecture, I give an introduction to one of the most active and exciting 
recent developments in theoretical particle physics: the discovery of heavy quark 
symmetry and the development of the heavy quark effective field theory. For a 
more detailed description, I refer to my recent review article in Physics Reports^. 
Earlier introductions into this field were given by Georgi 
WiseEil, and ManneE. 
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GrinsteinEIHI, Isgur and 



3.1. The Physical PicturM^ 

There are several reasons why the strong interactions of systems containing 
heavy quarks are easier to understand than those of systems containing only light 
quarks. The first is asymptotic freedom, the fact that the effective coupling constant 
of QCD becomes weak in processes with large momentum transfer, corresponding 
to interactions at short-distance scalesllil. At large distances, on the other hand, 
the coupling becomes strong, leading to nonperturbative phenomena such as the 
confinement of quarks and gluons on a length scale i?had ~ 1/Aqcd ~ 1 &n, which 
determines the typical size of hadrons. Roughly speaking, Aqcd ~ 0.2 GeV is the 
energy scale that separates the regions of large and small coupling constant. When 
the mass of a quark Q is much larger than this scale, ttiq ^> Aqcd, Q is called a heavy 
quark. The quarks of the standard model fall naturally into two classes: up, down 
and strange are light quarks, whereas charm, bottom and top are heavy quarks? 
For heavy quarks, the effective coupling constant Q s (mg) is small, implying that 
on length scales comparable to the Compton wavelength Xq ~ the strong 

interactions are perturbative and much like the electromagnetic interactions. In 
fact, the quarkonium systems (QQ), whose size is of order Xq/ a s {mo) <C -Rhad, 
are very much hydrogen-like. Since the discovery of asymptotic freedom, their 
properties could be predictedEli before the observation of charmonium, and later of 
bottomonium states. 

Things are more complicated for systems composed of a heavy quark and other 
light constituents. The size of such systems is determined by -Rhad, and the typical 
momenta exchanged between the heavy and light constituents are of order Aqcd- 
The heavy quark is surrounded by a most complicated, strongly interacting cloud 
of light quarks, antiquarks, and gluons. This cloud is sometimes referred to as the 
"brown muck" , a term invented by Isgur to emphasize that the properties of such 
systems cannot be calculated from first principles (at least not in a perturbative 
way). In this case, it is the fact that rriQ ^> Aqcd, or better Xq <C -Rhad, i-e. the 
fact that the Compton wavelength of the heavy quark is much smaller than the size 
of the hadron, which leads to simplifications. To resolve the quantum numbers of 
the heavy quark would require a hard probe. The soft gluons which couple to the 

■^Ironically, the top quark is of no relevance to my discussion here, since it is too heavy to form 
hadronic bound states before it decays. 
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"brown muck" can only resolve distances much larger than Xq. Therefore, the light 
degrees of freedom are blind to the flavour (mass) and spin orientation of the heavy 
quark. They only experience its colour field, which extends over large distances 
because of confinement. In the rest frame of the heavy quark, it is in fact only the 
electric colour field that is important; relativistic effects such as colour magnetism 
vanish as tuq — > oo. Since the heavy quark spin participates in interactions only 
through such relativistic effects, it decouples. That the heavy quark mass becomes 
irrelevant can be seen as follows: As niq — > oo, the heavy quark and the hadron that 
contains it have the same velocity. In the rest frame of the hadron, the heavy quark 
is at rest, too. The wave function of the "brown muck" follows from a solution of 
the field equations of QCD subject to the boundary condition of a static triplet 
source of colour at the location of the heavy quark. This boundary condition is 
independent of uiq, and so is the solution for the configuration of the light degrees 
of freedom. 

It follows that, in the niq — > oo limit, hadronic systems which differ only in the 
flavour or spin quantum numbers of the heavy quark have the same configuration of 
their light degrees of freedom. Although this observation still does not allow us to 
calculate what this configuration is, it provides relations between the properties of 
such particles as the heavy mesons B, D, B* and D*, or the heavy baryons A& and 
A c (to the extent that corrections to the infinite quark mass limit are small in these 
systems) . Isgur and Wise realized that these relations result from new symmetries of 
the effective strong interactions of heavy quarks at low energies^. The configuration 
of light degrees of freedom in a hadron containing a single heavy quark Q(v, s) with 
velocity v and spin s does not change if this quark is replaced by another heavy 
quark Q'(v,s') with different flavour or spin, but with the same velocity. Both 
heavy quarks lead to the same static colour field. It is not necessary that the heavy 
quarks Q and Q' have similar masses. What is important is that their masses are 
large compared to Aqcd- For Nh heavy quark flavours, there is thus an SU(2Ah) 
spin-flavour symmetry group. These symmetries are in close correspondence to 
familiar properties of atoms: The flavour symmetry is analogous to the fact that 
different isotopes have the same chemistry, since to good approximation the wave 
function of the electrons is independent of the mass of the nucleus. The electrons 
only see the total nuclear charge. The spin symmetry is analogous to the fact that 
the hyperfine levels in atoms are nearly degenerate. The nuclear spin decouples in 
the limit m e /rriN —* 0. 

The heavy quark symmetry is an approximate symmetry, and corrections arise 
since the quark masses are not infinite. These corrections are of order Aqcd/^q- 
The condition niq ^> Aqcd is necessary and sufficient for a system containing a 
heavy quark to be close to the symmetry limit. In many respects, heavy quark 
symmetry is complementary to chiral symmetry, which arises in the opposite limit 
of small quark masses, m q <C Aq CD . There is an important distinction, however. 
Whereas chiral symmetry is a symmetry of the QCD Lagrangian in the limit of 
vanishing quark masses, heavy quark symmetry is not a symmetry of the Lagrangian 
(not even an approximate one), but rather a symmetry of an effective theory which 
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is a good approximation to QCD in a certain kinematic region. It is realized only 
in systems in which a heavy quark interacts predominantly by the exchange of soft 
gluons. In such systems the heavy quark is almost on-shell; its momentum fluctuates 
around the mass shell by an amount of order Aqcd- The corresponding changes in 
the velocity of the heavy quark vanish as Aqcd/^q ~~ > 0. The velocity becomes a 
conserved quantity and is no longer a dynamical degree of freedomEiJ. 

Results derived based on heavy quark symmetry are model-independent conse- 
quences of QCD in a well-defined limit. The symmetry-breaking corrections can, 
at least in principle, be studied in a systematic way. To this end, it is however 
necessary to recast the QCD Lagrangian for a heavy quark, 

£ Q = Q(zp-m Q )Q, (99) 

into a form suitable for taking the limit tuq — > oo. 

3.2. Heavy Quark Effective TheoryE-E 

In particle physics, it is often the case that the effects of a very heavy particle 
become irrelevant at low energies. It is then useful to construct a low-energy effective 
theory, in which this heavy particle no longer appears. Eventually, this effective 
theory will be easier to deal with than the full theory. A familiar example is Fermi's 
theory of the weak interactions. For the description of weak decays of hadrons, 
one can approximate the weak interactions by point-like four-fermion couplings 
governed by a dimensionful coupling constant Gp. Only at energies much larger 
than the masses of hadrons can one resolve the structure of the intermediate vector 
bosons W and Z. 

Technically, the process_of removing the degrees of freedom of a heavy particle 
involves the following steps&Z^ One first identifies the heavy particle fields and 
"integrates them out" in the generating functional of the Green functions of the 
theory. This is possible since at low energies the heavy particle does not appear 
as an external state. However, although the action of the full theory is usually a 
local one, what results after this first step is a nonlocal effective action. The non- 
locality is related to the fact that in the full theory the heavy particle (with mass 
M) can appear in virtual processes and propagate over a short but finite distance 
Ax ~ 1/M. Thus a second step is required to get to a local effective Lagrangian: 
The nonlocal effective action is rewritten as an infinite series of local terms in an 
operator product expansion (OPE)Sa Roughly speaking, this corresponds to an 
expansion in powers of 1/M. It is in this step that the short- and long-distance 
physics is disentangled. The long-distance physics corresponds to interactions at 
low energies and is the same in the full and the effective theory. But short-distance 
effects arising from quantum corrections involving large virtual momenta (of order 
M) are not reproduced in the effective theory, once the heavy particle has been 
integrated out. In a third step, they have to be added in a perturbative way using 
renormalization-group techniques. This procedure is called "matching" . It leads to 
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renormalizations of the coefficients of the local operators in the effective Lagrangian. 
An example is the effective Lagrangian for nonleptonic weak decays, in which radia- 
tive corrections from hard gluons with virtual momenta in the range between my/ 
and some renormalization scale \i ~ f GeV g ive rise to Wilson coefficients, which 



renormalize the local four-fermion interactionstlHrt^ . 

The heavy quark effective theory (HQET) is constructed to provide a simplified 
description of processes where a heavy quark interacts with light degrees of freedom 
by the exchange of soft gluons. Clearly, tuq is the high-energy scale in this case, 
and Aqcd is the scale of the hadronic physics one is interested in. However, a 
subtlety arises since one wants to describe the properties and decays of hadrons 
which contain a heavy quark. Hence, it is not possible to remove the heavy quark 
completely from the effective theory. What is possible, however, is to integrate out 
the "small components" in the full heavy quark spinor, which describe fluctuations 
around the mass shell. 

The starting point in the construction of the low-energy effective theory is the 
observation that a very heavy quark bound inside a hadron moves with essentially 
the hadron's velocity v and is almost on-shell. Its momentum can be written as 

p% = m Q v» + k», (fOO) 

where the components of the so-called residual momentum /c M are much smaller 
than rriQ. Interactions of the heavy quark with light degrees of freedom change the 
residual momentum by an amount of order Ak^ ~ Aqcd, but the corresponding 
changes in the heavy quark velocity vanish as Aqcd/^q — > 0. In this situation, it 
is appropriate to introduce "large" and "small" component fields h v and H v by 

h v ( x ) = e im Q v - x P + Q(x), 

H v (x) =e im ^P_Q(x), (101) 
where P + and P_ are projection operators defined as 

P ± = hM. (10 2) 

It follows that 

Q(x) = e - im ^ x [h v (x) + H v (x)] . (103) 

Because of the projection operators, the new fields satisfy fh v = h v and fti v = —H v . 
In the rest frame, h v corresponds to the upper two components of Q, while H v 
corresponds to the lower ones. Whereas h v annihilates a heavy quark with velocity 
v, H v creates a heavy antiquark with velocity v. If the heavy quark was on-shell, 
the field H v would be absent. 

In terms of the new fields, the QCD Lagrangian (|99| ) for a heavy quark takes 
the following form: 

C Q = h v iv-Dh v - H v (iv-D + 2m Q ) H v + h v ift ± H v + H v ip±h v , (104) 
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where Dj_ = D^ 1 — t> M v-D is orthogonal to the heavy quark velocity: v-D± = 0. In the 
rest frame, D± = (0, -D ) contains the spatial components of the covariant derivative. 
From ( |104|) , it is apparent that h v describes massless degrees of freedom, whereas H v 
corresponds to fluctuations with twice the heavy quark mass. These are the heavy 
degrees of freedom, which will be eliminated in the construction of the effective 
theory. The fields are mixed by the presence of the third and fourth terms, which 
describe pair creation or annihilation of heavy quarks and antiquarks. As shown in 



the first diagram in Fig. [Fy, in a virtual process a heavy quark propagating forward 
in time can turn into a virtual antiquark propagating backward in time, and then 
turn back into a quark. The energy of the intermediate quantum state hhH is 
larger than the energy of the incoming heavy quark by at least 2mg. Because of 
this large energy gap, the virtual quantum fluctuation can only propagate over a 
short distance Ax ~ 1/^q- On hadronic scales set by Rh&d = 1/Aqcd, the process 
essentially looks like a local interaction of the form 

h v ip ±7 —ip x h v , (105) 

where I have simply replaced the propagator for H v by l/2rriQ. A more correct treat- 
ment is to integrate out the small component field H v , thereby deriving a nonlocal 
effective action for the large component field h v , which can then be expanded in 
terms of local operators. Before doing this, let me mention a second type of virtual 
corrections that involve pair creation, namely heavy quark loops. An example is 
depicted in the second diagram in Fig. O. Heavy quark loops cannot be described 
in terms of the effective fields h v and H v , since the quark velocities in the loop are 
not conserved, and are in no way related to hadron velocities. However, such short- 
distance processes are proportional to the small coupling constant a s (mg) and can 
be calculated in perturbation theory. They lead to corrections which are added onto 
the low-energy effective theory in the matching procedure. 




Figure 11: Virtual fluctuations involving pair creation of 
heavy quarks. Time flows to the right. 

On a classical level, the heavy degrees of freedom represented by H v can be 
eliminated using the equations of motion of QCD. Substituting (|103|) into {i p — 

m o) Q — gives 

iph v + {ifl-2m Q )H v = Q, (106) 
and multiplying this by P± one derives the two equations 

— iv-D h v = i p±H v , 
(iv-D + 2m Q )H v = ifl ± h v . (107) 
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The second can be solved to give (77 — > +0) 

Hv = To 7 1 „ • S ip±K, (108) 

(2mg + iv -D — iT]) 

which shows that the small component field H v is indeed of order I/tuq. One can 
now insert this solution into the first equation to obtain the equation of motion for 
h v . It is easy to see that this equation follows from the nonlocal effective Lagrangian 

C cS = h v iv-Dh v + h v i]p ± — 1 — ip_i_h v . (109) 

(2mQ + iv -D — irj) 

Clearly, the second term precisely corresponds to the first class of virtual processes 



depicted in Fig. |11 



Mannel, Roberts, and Ryzak have derived this Lagrangian in a more elegant way 
by manipulating the generating functional for QCD Green's functions containing 
heavy quark fields^. They start from the field redefinition ( |103|) and couple the 



large component fields h v to external sources p v . Green's functions with an arbitrary 
number of h v fields can be constructed by taking derivatives with respect to p v . 
No sources are needed for the heavy degrees of freedom represented by H v . The 
functional integral over these fields is Gaussian and can be performed explicitly, 
leading to the nonlocal effective action 

S cS = Jd 4 x£ eS -i\nA, (110) 

with C e s as given in (|109| ). The appearance of the logarithm of the determinant 

A = exp (— Tr In [2m,Q + iv-D — irj\\ (HI) 

is a quantum effect not present in the classical derivation given above. However, in 
this case the determinant can be regulated in a gauge-invariant way, and by choosing 
the axial gauge v ■ A = one shows that In A is just an irrelevant constantll^lli^. 

Because of the phase factor in ( |1U3| ), the x dependence of the effective heavy 
quark field h v is weak. In momentum space, derivatives acting on h v produce 
powers of the residual momentum k, which is much smaller than mq. Hence, the 
nonlocal effective Lagrangian (|1U9| ) allows for a derivative expansion in powers of 
iD/rriQ. Taking into account that h v contains a P + projection operator, and using 
the identity 



P + ip ± ilp ± P + = P+ 



(^x) 2 + |a a ,G^ 



P+ , (112) 



where [iD a ,iD@] = ig s G a/3 is the gluon field-strength tensor, one finds that@0 

£ eff = h v iv-Dh v + h v {iD ± ) 2 h v + h v a aP G af3 h v + 0(l/m Q ) . (113) 

Iraq ^rriQ 
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j = igCTa^U 01 



a, a 



Figure 12: Feynman rules of the heavy quark effective the- 
ory and a are colour indices). It has become standard 
to represent a heavy quark by a double line. The residual 
momentum k is defined in fllOOj ). 

In the limit tuq — > oo, only the first terms remains: 

Coo — h v iv-D h v . (114) 

This is the effective Lagrangian of HQET. From it are derived the Feynman rules 
depicted in Fig. IT2]. Let me take a moment to study the symmetries of this 



LagrangiaD^H Since there appear no Dirac matrices, interactions of the heavy quark 
with gluons leave its spin unchanged. Associated with this is an SU(2) symmetry 
group, under which is invariant. The action of this symmetry on the heavy 
quark fields becomes most transparent in the rest frame, where the generators S l 
of spin SU(2) can be chosen as 

where a 1 are the Pauli matrices. In a general frame, one defines a set of three 
orthonormal vectors e l orthogonal to v, and takes the generators of the spin sym- 
metry as S l = i % ■ These matrices satisfy the commutation relations of SU(2) 
and commute with ^: 

[S i ,S j ]=ie ijk S k , [^S l }=0. (116) 

An infinitesimal SU(2) transformation h v — > (1 + ie- S) h v leaves the Lagrangian 
invariant, 

6Coo — h v [iv-D, ie ■ S] h v — , (H7) 
and preserves the on-shell condition $h v = h v , since 

f{l+ie-§)h v = {l + ie- S ) j h v = (1 + ie ■ S ) h v . (118) 

There is another symmetry of HQET, which arises since the mass of the heavy 
quark does not appear in the effective Lagrangian. When there are Nh heavy quarks 
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moving at the same velocity, one can simply extend ( |114| ) by writing 

N h _ 

C oc = J2K^-DK. (119) 

1=1 

This is clearly invariant under rotations in flavour space. When combined with the 
spin symmetry, the symmetry grovm becomes promoted to SU(2iV/j). This is the 
heavy quark spin-flavour symmetryHaEil. Its physical content is that, in the vtlq — ► oo 
limit, the strong interactions of a heavy quark become independent of its mass and 
spin. 

Let me now have a look at the operators appearing at order I/trq in ( |113|) . 
They are easiest to identify in the rest frame. The first operator, 

Ckin = 7T—hv ( W ±) 2 K^-TT—K (iD f h v , (120) 
2m Q 2m Q 

is the gauge-covariant extension of the kinetic energy arising from the off-shell resid- 
ual motion of the heavy quark. The second operator is the non-Abelian analogue of 
the Pauli term, which describes the colour-magnetic interaction of the heavy quark 
spin with the gluon field: 

Cmag = ~~r~~ h v a aj3 G a/3 h v — > : — h v S-B c h v . (121) 

4m Q m Q 



Here S is the spin operator defined in (|115|) , and B l c = — |e ufc G jfe are the components 



of the colour-magnetic field. This chromo-magnetic interaction is a relativistic effect, 
which scales like 1 / iriq . This is the origin of the heavy quark spin symmetry. 

3.3. Spectroscopic Implications 

The spin-flavour symmetry leads to many interesting relations between the prop- 
erties of hadrons containing a heavy quark. The most direct consequences concern 



the spectroscopy of such states^. In the mg — > oo limit, the spin of the heavy 
quark and the total angular momentum j of the light degrees of freedom are sepa- 
rately conserved by the strong interactions. Because of heavy quark symmetry, the 
dynamics is independent of the spin and mass of the heavy quark. Hadronic states 
can thus be classified by the quantum numbers (flavour, spin, parity, etc.) of the 
light degrees of freedom. The spin symmetry predicts that, for fixed j ^ 0, there 
is a doublet of degenerate states with total spin J = j ± -. The flavour symmetry 
relates the properties of states with different heavy quark flavour. 

Consider, as an example, the ground-state mesons containing a heavy quark. In 
this case the light degrees of freedom have the quantum numbers of an antiquark, 
and the degenerate states are the pseudoscalar ( J = 0) and vector (J = 1) mesons. 
In the charm and bottom systems, one knows experimentally 

rriB* — rn B ~ 46 MeV , 

m D * - m D ~ 142 MeV . (122) 
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These mass splittings are in fact reasonably small. To be more specific, at order 
I/wig one expects hyperfine corrections to resolve the degeneracy, for instance tub* — 
itlb oc l/m&. This leads to the refined prediction m B * — m B ~ m 2 D , — m 2 D ~ const. 
The data are compatible with this: 

m|, - m| ~ 0.49 GeV 2 , 

m 2 D , -m 2 D ~ 0.55 GeV 2 . (123) 

One can also study excited meson states, in which the light constituents carry 
orbital angular momentum. It is tempting to interpret _Di(2420) with J p = 1 + and 
-D2(2460) with J p = 2 + as the spin doublet corresponding to j = |. The small mass 
difference mm — ~ 35 MeV supports this assertion. One then expects 



m l* ~ ™> 2 Bl ^ m 2 D . - m 2 Di ~ 0.17 GeV 2 (124) 



for the corresponding states in the bottom system. The fact that this mass splitting 
is smaller than for the ground-state mesons is not unexpected. For instance, in the 
nonrelativistic constituent quark model the light antiquark in these excited mesons 
is in a p-wave state, and its wave function at the location of the heavy quark 
vanishes. Hence, in this model hyperfine corrections are strongly suppressed. 

A typical prediction of the flavour symmetry is that the "excitation energies" 
for states with different quantum numbers of the light degrees of freedom are ap- 
proximately the same in the charm and bottom systems. For instance, one expects 

rriB s — tub — rri£) s — rn D ~ 100 MeV , 
rriB 1 — rriB — tud^ — rnD — 557 MeV , 

n^B* — n^B — ran* — mo — 593 MeV . (125) 

The first relation has been very nicely confirmed by the discovery of the B$ meson 
at The observed masalHJ m# s = 5.368 ± 0.005 GeV corresponds to an 



excitation energy m^- m^ = 89±5 MeV. 



3.4- Weak Decay Form Factors 

Of particular interest are the relations between the weak decay form factors of 
heavy mesons, which parametrize hadronic matrix elements of currents between two 
meson states containing a heavy quark. These relations have been derived by Isgur 
and WiseEl, generalizing ideas developed by Nussinov and WetzelS and by Voloshin 
and ShifmanE-M. For the purpose of this discussion, it is convenient to work with 
a mass-independent normalization of meson states, 

(M(p')\M(p)) = -E- (2n) 3 5 3 (p-p') , (126) 
m M 

instead of the more conventional, relativistic normalization 

(M(p')\M(p)} = 2p°(27r) 3 5 3 (p-p'). (127) 
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In the first case, p /rriM = v° depends only on the meson velocity. In fact, it is more 
natural for heavy quark systems to use velocity rather than momentum variables. 
I will thus write \M(v)) instead of \M(p)). The relation to the conventionally 

— 1/2 — ' 

normalized states is \M(v)) = m M M(p)). 



t < t t > t 

Figure 13: The action of an external heavy quark current, 
as seen by the light degrees of freedom in the initial state. 

Consider now the elastic scattering of a pseudoscalar meson, P{y) — > P(v'), 
induced by an external vector current coupled to the heavy quark contained in P. 
Before the action of the current, the light degrees of freedom in the initial state 
orbit around the heavy quark, which acts as a source of colour moving with the 
meson's velocity v. On average, this is also the velocity of the "brown muck". The 
action of the current is to replace instantaneously (at t = to) the colour source by 
one moving at velocity v ', as indicated in Fig. [L3[ If v — v', nothing really happens; 
the light degrees of freedom do not realize that there was a current acting on the 
heavy quark. If the velocities are different, however, the "brown muck" suddenly 
finds itself interacting with a moving (relative to its rest frame) colour source. Soft 
gluons have to be exchanged in order to rearrange the light degrees of freedom 
and build the final state meson moving at velocity v' . This rearrangement leads 
to a form factor suppression, which reflects the fact that as the velocities become 
more and more different, the probability for an elastic transition decreases. The 
important observation is that, in the tuq — > oo limit, the form factor can only 
depend on the Lorentz boost 7 = v ■ v' that connects the rest frames of the initial- 
and final-state mesons. That the form factor is independent of the heavy quark 
mass can also be seen from the following intuitive argument: The light constituents 
in the initial and final state carry momenta that are typically of order Aqcd^ and 
Aqcd^', respectively. Thus, their momentum transfer is q 2 ~ AQ CD (t> • v' — 1), which 
is in fact independent of tuq. 

The result of this discussion is that in the effective theory, which provides the 
appropriate framework to consider the limit uiq — ► 00 with the quark velocities kept 
fixed, the hadronic matrix element describing the scattering process can be written 
asEl 

(P(v')\ h v ,^h v \P{v)) = i{v ■ v 1 ) (v + v'Y , (128) 
with a form factor £(v ■ v') that does not depend on tuq. Since the matrix element 
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is invariant under complex conjugation combined with an interchange of v and t/, 
the function £(t> • v') must be real. That there is no term proportional to (v — v')^ 
can be seen by contracting the matrix element with (v — t>') M , and using i/ih v = h v 
and h v if = h v /. 

One can now use the flavour symmetry to replace the heavy quark Q in one 
of the meson states by a heavy quark Q' of a different flavour, thereby turning 
P into another pseudoscalar meson P'. At the same time, the current becomes 
a flavour-changing vector current. In the infinite mass limit, this is a symmetry 
transformation, under which the effective Lagrangian is invariant. Hence, the matrix 
element 

(P V) I Kl^K \P(v)) = £(v • v') (v + v'Y (129) 

is still determined by the same function £(v -v'). This universal form factor is called 
the Isgur-Wise function, after the discoverers of this relationll 3 ! 

For equal velocities, the vector current J 71 = h' v ^h v = h' v v^h v is conserved in 
the effective theory, irrespective of the flavour of the heavy quarks, since 

id^Ji* = h' v (iv-Dh v ) + (iv-Dh' v ) h v = (130) 

by the equation of motion, iv-D h v = 0, which follows from the effective Lagrangian 
( |114j ). The conserved charges 

N Q/Q = J d 3 x J°(x) = J d 3 x h$(x) h v (x) (131) 

are generators of the flavour symmetry. The diagonal generators simply count the 
number of heavy quarks, whereas the off-diagonal ones replace a heavy quark by 
another: N QQ \P(v)) = \P{v)) and N QV \P{v)) = \P'(v)). It follows that 

(P'(v)\N Q , Q \P(v)) = (P(v)\P(v)) = 2v° (2tt) 3 5 3 (0 ) , (132) 

and from a comparison with (|129| ) one concludes that the Isgur-Wise function is 
normalized at the point of equal velocities: 

£(!) = !■ (133) 

This can easily be understood in terms of the physical picture discussed above: 
When there is no velocity change, the light degrees of freedom see the same colour 
field and are in an identical configuration before and after the action of the current. 
There is no form factor suppression. Since -E reC oii = m P' (v ■ v' — 1) is the recoil 
energy of the daughter meson P' in the rest frame of the parent meson P, the point 
v ■ v' = 1 is referred to as the zero recoil limit. 

The heavy quark spin symmetry leads to additional relations among weak decay 
form factors. It can be used to relate matrix elements involving vector mesons to 
those involving pseudoscalar mesons. In the effective theory, a vector meson with 
longitudinal polarization e 3 is related to a pseudoscalar meson by 

\V(v,e 3 )) = 2S 3 Q \P(v)), (134) 
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where Sq is a Hermitian operator that acts on the spin of the heavy quark Q. It 
follows that 



(V'(v',e 3 )\K,Th v \P(v)) = (P'(v')\2[S 3 Qi rK^h v )\P(v)) 

= (P'(v')\h' v ,(2S 3 T)h v \P(v)), (135) 

where T can be an arbitrary combination of Dirac matrices, and S 3 is a matrix 
representation of the operator Sq,. It is easiest to evaluate this expression in the 
rest frame of the final-state meson, where [cf. (|115| )1 

^ = (1,0,0,0), e£ = (0,0,0,l), 5 3 = ^ 7 57°7 3 . (136) 

It is then straightforward to obtain the following commutation relations for the 
components of the weak current (V — A)^ = /i^,7 M (l — 75) h v : 

2 [S Q ,, V° - A } = A 3 -V 3 , 
2 [Sq,, V 3 - A 3 ] = A°-V°, 
2[S Q/ ,V 1 -A 1 ]= t(A 2 -V 2 ), 

2 [Sq,, V 2 — A 2 ] = —i^A 1 — V 1 ) . (137) 

Using ( |135| ) and ( |137| ), one can relate the matrix element of the weak current be- 
tween a pseudoscalar and a vector meson to the matrix element of the vector current 
between two pseudoscalar mesons given in ( |129| ). The result can be written in the 
covariant formE^I: 

(V'(v',e)\ h' v ,r(l - 75) K \P(v)) = k^e» e(« • «0 

^ {vv' + l)- v'" e* ■ v] £(y ■ v') , (138) 

where e 0123 = —£0123 = — 1- Once again, the matrix element is completely described 
in terms of the universal Isgur-Wise form factor. 

Equations (|129| ) and (|138| ) summarize the relations imposed by heavy quark sym- 
metry on the weak decay form factors describing the semileptonic decay processes 
B — > D Iv and B — > D*i v. These relations are model-independent consequences of 
QCD in the limit where rrib,m c ^> Aqcd- They play a crucial role in the determi- 
nation of I Vrf,|, as I will discuss later in this lecture. 



3.5. The 1/rriQ Expansion 

At tree level, eq. ( p.!3| ) defines the Lagrangian of HQET as a series of local, 
higher-dimension operators multiplied by powers of 1/toq. The expression for H v 
in (|108| ) can be used to derive a similar expansion for the full heavy quark field 
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Q(x) 



Q(x) = e 



ITTIQV-X 



1 + 7 ^ \ i Pa 

{2rriQ + iv-D — irj) 

ip 



h v (x) 



= e- im ^ x (l + i^i + ...)h v (x). (139) 
V 2m Q J 

This relation contains the recipe for the construction (at tree level) of any operator 
in HQET that contains one or more heavy quark fields. For instance, the vector 
current V 1 = q'j^Q composed of a heavy and a light quark is represented as 

V»{x) = e~ im ^ x q(x) Y ( 1 + ^ + . . . ) h v {x) . (140) 

V 2m Q ) 

Matrix elements of this current can be parametrized by hadronic form factors, and 
the purpose of using an effective theory is to make the tuq dependence of these form 
factors explicit. 

Consider, as an example, the matrix element of V^(0) between a heavy meson 
M(v) and some light final state £: 

(£\V»\M(v)) = (e\q-fh v \M(v)) + -—(£\q'fip ± h v \M(v)) + ... . (141) 

2rriQ 

It would be nice if the matrix elements on the right-hand side of this equation were 
independent of tjiq. Then the second term would give the leading power correction 
to the first one. However, the equation of motion for h v derived from ( |113|) contains 
1/rriQ corrections, too. This leads to an itlq dependence of any hadronic matrix 
element of operators containing such fields. Another way to say this is that the 
eigenstates of the effective Lagrangian C cS (supplemented by the standard QCD 
Lagrangian for the light quarks and gluons) depend, at higher order in 1/toq, on the 
heavy quark mass. This is no surprise, since the effective Lagrangian is equivalent 
to the Lagrangian of the full theory. 

It is better to organize things in a slightly different way by working with the 
eigenstates of only the leading-order effective Lagrangian in ( 11 14T) . t reating the 
higher-dimension operators perturbatively as power correctiongBS 184 !! 187 !. Then the 
equation of motion satisfied by h v is simply 

iv-Dh v = 0, (142) 

and the states of the effective theory are truly independent of tuq. However, these 
states are now different from the states of the full theory. For instance, instead of 
(11411) one should write 



+ 7 ^—(e\qi fl tp±h v \M(v)) c 



(£\V>\M(v)) QpD = (£\q>fh v \M(v)) 

1 

2m. 
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' ( i 1 i I dy T { q 7^(0), h v (iD L ) 2 h v (y) } \M(v)) c 



2m^ 

+ (I | i [dy T { q^h v (0), h v a aP G a ? h v (y) } \M(v))oq 

+ 0(l/m 2 Q ). (143) 

The matrix elements on the right-hand side are independent of the heavy quark 
mass. The mass dependence of the states \M(v))qcd is reflected in HQET by the 
appearance of the third and fourth term, which arise from an insertion of the first- 
order power corrections in the Lagrangian into the leading-order matrix element 
of the current. This insertion can be thought of as being a correction to the wave 
function of the heavy meson. From now on I will omit the subscripts on the states. 
I will instead use the symbol "=" and write 

^(0) = g 7 ^(0) + J—q^ip x h v (0) 

+ 7^— ifdyT{ q^h v (0), h v {iD ± ) 2 h v {y) } 
2m Q J 

+ -^-i fdy T {q^h v (0),h v v a pG al3 h v (y)} + . . . (144) 
4m<2 J 

to indicate that the operators on the two sides of this equation have to be evaluated 
between different states. 

3. 6. Matching and Running 

In Section 3.2, I have discussed the first two steps in the construction of HQET. 
Integrating out the small components in the heavy quark spinor fields, a nonlocal 
effective action was derived. This was then expanded in a series of local operators of 
increasing dimension to obtain an effective Lagrangian. A similar expansion could 
be written down for any external current. The effective Lagrangian and the effective 
currents derived that way correctly reproduce the long-distance physics of the full 
theory. They cannot describe the short-distance physics correctly, however. The 
reason is obvious: The heavy quark participates in strong interactions through its 
coupling to gluons. These gluons can be soft or hard, i.e. their virtual momenta can 
be small, of the order of the confinement scale, or large, of the order of the heavy 
quark mass. But hard gluons can resolve the nonlocality of the propagator of the 
small component fields H v . Their effects are not taken into account in the naive 
version of the OPE, which was used in the derivation of the effective Lagrangian in 
(|113| ) and the effective vector current in ( |144| ). So far, the effective theory provides 



an appropriate description only at scales /i ttiq. 

In this section, I will discuss the systematic treatment of short-distance cor- 
rections. A new feature of such corrections is that through the running coupling 
constant they induce a logarithmic dependence on the heavy quark mass0. The 
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important observation is that a s (mo) is small, so that these effects can be cal- 
culated in perturbation theory. Consider, as an example, matrix elements of the 
vector current V* = O'y^Q- In QCD this current is partially conserved and needs 
no renormalizationE^l its matrix elements are free of ultraviolet divergences. Still, 
these matrix elements can have logarithmic dependence on tuq from the exchange 
of hard gluons with virtual momenta of the order of the heavy quark mass. If one 
goes over to the effective theory by taking the limit tuq — > oo, these logarithms 
diverge. Consequently, the vector current in the effective theory does require a 
renormalizationlia. Its matrix elements depend on an arbitrary renormalization 
scale /x, which separates the regions of short- and long-distance physics. If /i is 
chosen such that Aq CD <C /i <C mg, the effective coupling constant in the region 
between /i and tuq is small, and perturbation theory can be used to compute the 
short-distance corrections. These corrections have to be added to the matrix ele- 
ments of the effective theory, which only contain the long-distance physics below 
the scale fi. Schematically, then, the relation between matrix elements in the full 
and in the effective theory is 

( VHmq) >qcd = C (m Q , /i) (V (aO>oo + g^p^ (V^))^ + . . . , (145) 

Q 

where I have indicated that matrix elements of V 11 in the full theory depend on mq , 
whereas matrix elements of operators in the effective theory are mass-independent, 
but do depend on the renormalization scale. The Wilson coefficients Ci(m,Q,fj,) 
are defined by this relation. Order by order in perturbation theory, they can be 
computed from a comparison of the matrix elements in both theories. Since the 
effective theory is constructed to reproduce correctly the low-energy behaviour of the 
full theory, this "matching" procedure is independent of any long-distance physics, 
such as infrared singularities, nonperturbative effects, the nature of the external 
states used in the matrix elements, physical cuts, etc. Only at high energies do 
the two theories differ, and these differences are corrected for by the short-distance 
coefficients. 

The calculation of the coefficient functions in perturbation theory uses the pow- 
erful methods of the renormalization group, which I have discussed in my first 
lecture. It is in principle straightforward, yet in practice rather tedious. Much of 
the recent work on heavy quark symmetry has been devoted to this subject. In 
this section, I will discuss as an illustration the wave-function renormalization of 
the heavy quark field in HQET, and the short-distance expansion of a heavy-light 
vector current to leading order in l/m,Q. For a more comprehensive presentation of 
short-distance corrections, I refer to my review articleEl 

In quantum field theory, the parameters and fields of the Lagrangian have no 
direct physical significance. They have to be renormalized before they can be 
related to observable quantities. In an intermediate step the theory has to be 
regularized. The most convenient regularization scheme in QCD is dimensional 
regularizationHlHl, in which the dimension of space-time is analytically continued 
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to D = 4 — 2e, with e being infinitesimal. Loop integrals that are logarithmically 
divergent in four dimensions become finite for e > 0. From the fact that the action 
S = J d D x£(x) is dimensionless, one can derive the mass dimensions of the fields 
and parameters of the theory. For instance, one finds that the "bare" coupling con- 
stant a!? are is no longer dimensionless if D ^ 4: dim[a^ are ] = 2e. In a renormalizable 
theory, it is possible to rewrite the Lagrangian in terms of renormalized quantities, 
in such a way that Green's functions of the renormalized fields remain finite as 

, 1/2 

e — > 0. For QCD, one introduces renormalized quantities by Q e = ZA Q ven , 

A hare = Z]/ 2 A ven , a^ are = p? e Z a a v s en , etc., where /i is an arbitrary mass scale in- 
troduced to render the renormalized coupling constant dimensionless. Similarly, in 
HQET one defines a renormalized heavy quark field by h^ aie = Z]j 2 h T ^ n . From now 
on, the superscript "ren" will be omitted. 

As a warm-up, let me sketch the calculation of the wave-function renormalization 
for a heavy quark in HQET. In the MS scheme, Zh can be computed from the 1/e- 
pole in the quark self-energy shown in Fig. [TJ|: 

As long as v ■ k < 0, the self-energy is infrared finite and real. The result is gauge- 
dependent, however. In Feynman gauge, one obtains at one- loop order 



E{v ■ k) 



tigl r d D t 



(2vr) D {t 2 + i7])[v [t + k)+ir]} 



f dp r a D t 



3 J r J (2n) D { t 2 + 2pv(t + k)+ir ] } 2 

^ 

where p is a dimensionful Feynman parameter, and to = —2v ■ k > acts as an 
infrared cutoff. A straightforward calculation leads to 



4a, 



r(2e)r(l-e) (-^) 
3tt V ' V ; \ATip?) 



148) 



where I have substituted p — u (1 — z)/z. From an expansion around e = 0, one 
obtains 

Z h = l + ^. (149) 
one 

This result was first derived by Politzer and WiseS. In the meantime, the calcula- 
tion was also done at the two-loop orderEHI — EH . 
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Figure 14: Self-energy — iH(v ■ k) of a heavy quark in HQET. 
The velocity v is conserved by the strong interactions. 

Similar to the fields and coupling constants, any composite operator built from 
quark and gluon fields may require a renormalization beyond that of its component 
fields. Such operators can be divided into three classes: gauge-invariant operators 
that do not vanish by the equations of motion (class-I), gauge-invariant operators 
that vanish by the equations of motion (class-II), and operators which are not gauge- 
invariant (class-Ill). In general, operators with the same dimension and quantum 
numbers mix under renormali zati on However, things simplify if one works with 



the background field techniqueli^ - ^, which serves for quantizing gauge theories, 
preserving explicit gauge invariance. This offers the advantage that a class-I opera- 
tor cannot mix with class-Ill operators, so that only gauge-invariant operators need 
be considered^ . Furthermore, class-II operators are irrelevant since their matrix 
elements vanish by the equations of motion. It it thus sufficient to consider class-I 
operators only. For a set {Oi} of n class-I operators that mix under renormaliza- 
tion, one defines aniixw matrix of renormalization factors by (summation over 
j is understood) C^ are = ZijOj, such that the matrix elements of the renormal- 
ized operators Oj remain finite as e —>■ 0. In contrast to the bare operators, the 
renormalized ones depend on the subtraction scale via the \i dependence of Z^: 

H^- O t = (li^-ZrA 0^ ve = -l ik O k , (150) 



d/LA 1 V~ dfi %3 ) 3 
where 

lik = Z^ 1 /i— Z jk (151) 

are called the anomalous dimensions. It is sometimes convenient to introduce a com- 
pact matrix notation, in which O is the vector of renormalized operators, Z is the 
matrix of renormalization factors, and 7 denotes the anomalous dimension matrix. 
Then the running of the renormalized operators is controlled by the renormalization- 
group equation (RGE) 

^ + 7)0 = 0. (152) 

As an example, consider the short-distance expansion of the heavy-light vector 
current = qj^Q to leading order in I/itlq. In the full theory, the vector cur- 
rent is not renormalized. But the effective current operators in HQET do require 
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renormalization. The general form of the short-distance expansion reads 



V»{m Q ) = Ci{m Q , fj) Oi(n) + 0(l/m Q ) 

= Ci(rriQ,n) Z^fr) O)™ + 0(l/m Q ) , (153) 

where I have indicated the \x dependence of the renormalized operators. This is the 
correct generalization of (|145| ) in the case of operator mixing. In general, a complete 



set of operators with the same quantum numbers appears on the right-hand side. 
In the case at hand, there are two such operators, namely 

O x = q^h v1 2 = qv»h v . (154) 

From ( |152p and the fact that the product Cj(m<g,/i) 0j(/x) must be //-independent, 
one can derive the RGE satisfied by the coefficient functions. It reads 

/^-f)(7(m Q ,//) = 0, (155) 

where 7* denotes the transposed anomalous dimension matrix, and I have collected 
the coefficients into a vector. The solution of the RGE proceeds as I described in 
my first lecture; however, it is technically more complicated in the case when there 
is operator mixing. One obtains 

C{m Q ,n) = U(/j,,m Q )C(mQ,m Q ) , (156) 

with the evolution matrix0@'0 

U(fi, m Q ) = T a exp J da . (157) 

Here "T a " means an ordering in the coupling constant such that the couplings in- 
crease from right to left (for fi < tuq). This is necessary since, in general, the anoma- 
lous dimension matrices at different values of a do not commute: [7(0:1), 7(0:2)] 7^ 0. 
In the present case, however, these complications are absent, since the anomalous 
dimension matrix turns out to be proportional to the unit matrix. 



I will now illustrate the solution of ( |156|) to next-to-leading order in renormaliza- 



tion-group improved perturbation theory. There are two ingredients to this calcu- 
lation. First, the Wilson coefficients at the high-energy scale /x = mq must be 
calculated to one-loop order, 

C(m Q , m Q ) = (7 (0 ) + C {1) + . . . , (158) 

by matching QCD onto the effective theory. Then, in order to control the running of 
the HQET operators, the anomalous dimension matrix 7 must be calculated to two- 
loop accuracy. Let me start with a comparison of the one-loop matrix elements of the 
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vector current in the full and in the effective theory. It is legitimate to perform the 
matching calculation with on-shell quark states. Then the matrix elements can be 
written as (V^) = u q T M uq, where the heavy quark spinor satisfies fiiQ — uq. There 
will be infrared divergences in the calculation, which I regulate by introducing a 
fictitious gluon mass A. At tree level, the vertex function in both theories is simply 
p/x _ The one-particle irreducible one-loop diagrams are shown in Fig. [15|. 
They have to be supplemented by an on-shell wave-function renormalization of the 
external lines. The complete one-loop vertex function in QCD isE3 

a, /, m 2 11 M „ 2a, 



As required by the nonrenormalization theorem for partially conserved currents, this 
result is gauge-independent and ultraviolet finite-iH. In the limit of degenerate quark 
masses, the vector is conserved. The space integral over its time component is the 
generator of a flavour symmetry. It cannot be renormalized. When the symmetry is 
softly broken by the presence of mass splittings, this is only relevant for small loop 
momenta but does not affect the ultraviolet region. Thus, there can only be a finite 
renormalization. Notice, however, that (|159| ) does contain an infrared singularity, 
because the matrix element was calculated using unphysical states. 





m Q v 



k = 



Figure 15: One- loop corrections to the matrix elements of 
the vector current in QCD and in HQET. The wavy line rep- 
resents the current. The external momenta are on-shell. 

In the effective theory, the one-loop vertex function for any heavy-light current 
q T h v is found to bec3 

+S(i+'4n)} r - « 

where the Dirac structure of T is completely arbitrary. It is generally true that 
these operators are renormalized multiplicatively and irrespectively of their Dirac 
structured! . Since there is no approximate flavour symmetry relating light and 
heavy quarks in the effective theory, it is not unexpected that the matrix elements 
of the bare currents are ultraviolet divergent. In the MS scheme, the currents are 
renormalized by a diagonal matrix Z with components 



a* 



,22 



2vre 



'12 



'21 







(161) 
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From ( |153[ ), it then follows that the renormalized one-loop vertex function is 



1 HQET 



2tt 



A 2 



6 



Ci(mQ,A*)7 /1 + C 2 (m Q ,fj) 



(162) 



Notice that the A dependence is the same as in ( |159|) . The short-distance coefficients, 
which follow from a comparison of the two results, are independent of the infrared 
regulator: 



Ci(m Q ,/i) 
C 2 (m Q ,/j,) 



1 + ^ 

7T 

2a« 



In 



3tt 



(163) 



In particular, one obtains for the matching conditions at /i 

Aot s {m Q ) 



Ci(m Q ,m Q ) 
C 2 (m Q ,m Q ) 



1 

3tt 

2a s {m Q ) 
3tt 



(164) 



This completes the first step of the calculation. 

The next step is to calculate the one- and two-loop coefficients of the anomalous 
dimension, as defined in (|i~3|). From ( |160| ), it follows that the anomalous dimension 
matrix is proportional to the unit matrix, with the one-loop coefficient 



7o 



-4. 



(165) 



This is the so-called hybrid anomalous dimension of heavy-light currents derived by 
Voloshin and ShifmanG-I In the context of the effective theory, it has been calculated 
by Polilzer and Wise@. The two-loop coefficient 71 has been calculated by Ji and 



MusolfE-^l, and by Broadhurst and GrozmE^. They obtain 



Ti 



254 56 



9 



27 



7T 



20 

IT 



(166) 



The final result for the Wilson coefficients can be derived using (|20|). It is possible 
to write the answer in the factorized form Ci(rriQ,n) = Ci(rriQ) K^fi), where the 
dependence on the two scales has been completely separated. One finds 



Ci(m Q ) = a s (m Q ) 



-2/A) 



1 + 



a 



i m Q) 



7T 



C 2 {m Q ) = [a s (m Q ) 



2a s (m Q ) 
3tt 



S 



71 



(167) 
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where 

o _ 7iA) ~ A7o _ o 153 - 1971/ 381 + 28tt 2 - 30^ 

8/5 2 (33 - 2n f ) 2 36 (33 - 2n f ) ' 1 J 

5.7. Flavour- Changing Heavy Quark Currents 

The OPE of currents in the effective theory becomes considerably more compli- 
cated in the case of flavour-changing currents. The weak current C7 M (1 — 75) b is of 
this form, however, so it is necessary to consider this case in detail. The complica- 
tions arise from the fact that there are two different heavy quark masses, roj and 
m c . Thus the calculation of the Wilson coefficients becomes a two-scale problem. 
In addition, the coefficients will depend in a nontrivial way on the velocity transfer 

W = V ■ V '. 

There are two obvious ways of performing the transition from QCD to an effec- 
tive low-energy theory in which both the bottom and the charm quarks are treated 
as heavy quarks (in the sense of HQET): The transition can either be done in a 
single step, or by considering first an intermediate theory with a static bottom 
quark, but a dynamical charm quark. The latter becomes heavy in a second step. 
If one could solve perturbation theory to all orders, both treatments would lead to 
the same results for the Wilson coefficients. The calculation differs in both cases, 
however, and the results also differ if the perturbation series is truncated. 

To see what the differences are, suppose first the two heavy quarks have similar 
masses, i.e. m& ~ m c ~ m, with m being some average mass. It is then natural to 
remove at the same time the dynamical degrees of freedom of both heavy quarks. 
Let me consider the cases of the vector current = c^b and of the axial vector 
current = 07^756 explicitly and denote them collectively by J M . Each of these 
currents obeys an expansion of the form 

3 

J M (m 6 ,m c ) = ^Cjim^m^m^w^) Jj(fi) + 0(1 /m) , (169) 

3=1 

where Jj = h c v , Tjh b v are local current operators in the effective theory with matrices 

ri = y, r 2 = u", r 3 = ^ (170) 

for the vector current, and 

ri= 7 M 75, r 2 = ^ 75 , r 3 = ^ 75 (m) 

for the axial vector current. In the above expression, I have indicated that matrix 
elements of the original current J M between states of the full theory depend on 
77ift and m c , whereas matrix elements of the operators Jj between states of the 
effective theory are mass-independent, but do depend on the renormalization scale 
ji. The short-distance coefficients are functions of the heavy quark masses, the 
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renormalization scale, and the matching scale m. The dependence on m would 
disappear if one could sum the perturbative series to all orders. The advantage of 
this first approach is its simplicity. At leading order in the 1/m expansion, only 
three current operators contribute. Matrix elements of higher-dimension operators 
are suppressed by powers of AqcdA™- The short-distance coefficients contain the 
dependence on mj and m c correctly to a given order in a s , via matching at = m. 
The only disadvantage is the residual dependence on the matching scale m, which 
arises when one calculates to finite order in perturbation theory. Although in a 
next-to-leading order calculation the scale in the leading anomalous scaling factor 
is determined, one has no control over the scale in the next-to-leading corrections 
proportional to a s (m). Since m b > m > m c , this introduces an uncertainty of order 
a 2 s ln(m 6 /m c ). 

The alternative approach is to consider first, as an intermediate step for m;, > 
fi > m c , an effective theory with only a heavy bottom quark. Denoting the effective 
current operators of dimension (3 + k) in this theory by J> {m C} /i), indicating that 
their matrix elements depend both on the charm quark mass and the renormaliza- 
tion scale, I write the short-distance expansion as an expansion in l/m^: 



OO -I 

1 — „\ % k )t 



J»(m b ,m c ) = E — E^i K^fK,^). (172) 

fc=0 m b i 

Since the velocity of the charm quark is still a dynamical degree of freedom, in this 
intermediate effective theory there are only two dimension-three operators, namely 

jf = ciy^ jf ] = cT 2 h b v . (173) 

A major complication arises from the fact that matrix elements of the higher- 
dimension operators J- k ^ will, in general, scale like m k c . This compensates the 
prefactor \jm\. Consequently, these operators cannot be neglected even at leading 
order in the heavy quark expansion. One would thus have to deal with an infinite 
number of operators in order to keep track of the full dependence on the heavy 
quark masses. Ignoring this difficulty for the moment, one may use ( |172| ) to scale 
the currents from /i = m;, down to = m c , where the matching is done onto the final 
effective theory with two heavy quarks. In this step, the operator basis collapses 
considerably. When terms of order AqcdA^q (here I use mg generically for m c or 
rrib) are neglected on the level of matrix elements, only the three operators Jj{n) in 
( |169| ) remain. Each of the operators of the intermediate theory has an expansion in 
terms of these operators, with coefficients Eiji 

3 

j\ k \m C) m c ) = Y, m c {E§\m c> w,ti) J 5 {p) + 0(l/m c )} . (174) 

3=1 

Combining this with ( |172| ), one obtains 



3 

J M (m b ,m c ) = ^2Cj(m b ,m c ,w,fi) Jj(fi) + 0{l/m Q ) , (175) 

3=1 
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with evolution coefficients: 



197 



x> /m r \k. 



C j (m b ,m c ,w,fj) = ^2(—j J2 D i ( m i» m c) 4 K 1 "'/ 1 ) • ( 176 ) 

k=0 mb i 

In this expression, the matching scale m of the first approach does not appear. The 
coefficients depend either on a s (m b ) or a s (m c ), i.e. the scaling in the intermediate 
region m b > /i > m c is properly taken into account. To achieve this, however, it 
would be necessary to consider an infinite number of operators in the intermediate 
effective theory. This is, of course, not manageable. It is important to realize, 
however, that the short-distance coefficients in ( |169|) and (|175|) must agree, and 



that this equality must hold order by order in an expansion in the mass ratio 
m c /m b . Using this fact, one can combine the two approaches into a consistent 
next-to-leading order calculation^. For details of the very tedious calculation, the 
reader is referred to the literaturellS^i^. Below, I will only give numerical results. 



The short-distance coefficients can be written in the factorized form 

CP(m b ,m c ,w,iJ,) = Cf\m b ,m c ,w) K hh (w , fi) . (177) 

I denote the coefficients for the vector current by and those for the axial vector 
current by Cf. The /^-dependent function K^(w,fi) is universal and normalized to 
unity at zero recoil: K^l,^) = 1. The renormalization-group invariant coefficients 

C 4 contain all dependence on the heavy quark masses. In Table [y, I give the 
numerical values of these coefficients, which are obtained using m b = 4.80 GeV and 
m c = 1.45 GeV for the heavy quark masses, as well as A^g = 0.25 GeV (for rif = 4) 
in the two-loop expression for a s (/i). The corresponding coupling constants are 
a s (m b ) ~ 0.20 and a s (m c ) ~ 0.32. 



Table 1: 


Short-distance 


coefficient 


s for b - 


-> c transitions. 


w 


Ci 


c 2 


c 3 


ci 


ci 


ci 


1.0 


1.136 


-0.085 


-0.021 


0.985 


-0.122 


0.042 


1.1 


1.107 


-0.080 


-0.021 


0.965 


-0.115 


0.040 


1.2 


1.081 


-0.077 


-0.020 


0.946 


-0.109 


0.038 


1.3 


1.056 


-0.073 


-0.019 


0.927 


-0.103 


0.036 


1.4 


1.033 


-0.070 


-0.018 


0.910 


-0.098 


0.035 


1.5 


1.011 


-0.067 


-0.018 


0.894 


-0.094 


0.033 


1.6 


0.991 


-0.064 


-0.017 


0.878 


-0.089 


0.032 


1.7 


0.972 


-0.062 


-0.017 


0.864 


-0.086 


0.031 


1.8 


0.953 


-0.059 


-0.016 


0.850 


-0.082 


0.030 



Of particular interest for the model-independent determination of | V cb \ to be 
discussed below is the value of the short-distance coefficient Cf(w) at zero recoil. 
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One finds 
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VA = Cl{l) = x^ \ 1 + 1.561 ^(m c )-a s (m b ) _ c^rnj /8 + ^ 

ix ix V3 1 — z 

'25 14 _ a/OK 1 



where x = a s (m c )/a s (rrib), and z = m c /m b . For A^jg = 0.25 ± 0.05 GeV and 
z = 0.30 ± 0.05, this yields ij A = 0.986 ± 0.006. 

Let me stop, at this point, the discussion of short-distance effects. There have 
been important calculations that 1 have no time to mention here, such as the renor- 
malization of the higher-dimension operators in the effective Lagrangia and 



the short-distance expansion of currents to next-to- leading order inc23 I /nig. A dis- 
cussion of these calculations can be found in my review article!!. 

3. 8. Covariant Representation of States 

The purpose of the OPE discussed in the previous sections was to disentangle the 
short-distance physics related to length scales set by the Compton wavelengths of the 
heavy quarks from confinement effects relevant at large distances. This procedure 
makes explicit the mq dependence of any Green's function of the full theory, which 
contains one or more heavy quark fields. Hadronic matrix elements of heavy quark 
currents have a 1/mq expansion as shown in (|145|) . The HQET matrix elements 
on the right-hand side of this equation contain the long-distance physics associated 
with the interactions of the cloud of light degrees of freedom among themselves and 
with the background colour-field provided by the heavy quarks. These hadronic 
quantities depend in a most complicated way on the "brown muck" quantum num- 
bers of the external states, the quantum numbers of the current, and on the heavy 
quark velocities. They are related to matrix elements of effective current operators 
in HQET. Recall that these matrix elements are independent of the heavy quark 
masses, if the states of the effective theory are taken to be the eigenstates of the 
leading-order Lagrangian in ( |114j ). Matrix elements evaluated using these states 



have a well-defined behaviour under spin-flavour symmetry transformations. When 
combined with the requirements of Lorentz covariance, restrictive constraints on 
their structure can be derived. I will now introduce a very elegant formalism due to 
Bjorken^i and Falk et al-S^il, which allows one to derive the general form of such 



matrix elements in a straightforward manner. The clue is to work with a covariant 
tensor representation of states with definite transformation properties under the 
Lorentz group and the heavy quark spin-flavour symmetry. 

The eigenstates of HQET can be thought of as the "would-be hadrons" built 
from an infinitely heavy quark dressed with light quarks, antiquarks and gluons. In 
such a state, both the heavy quark and the cloud of light degrees of freedom have 
well-defined transformation properties under the Lorentz group. The heavy quark 
can be represented by a spinor Uh(v,s) satisfying ^JUhiy.s) = Uh(v,s), where v is 
the velocity of the hadron. Because of heavy quark symmetry, the wave function 
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of the state (when properly normalized) is independent of the flavour and spin of 
the heavy quark, and the states can be characterized by the quantum numbers of 
the "brown muck" . In particular, for each configuration of light degrees of freedom 
with total angular momentum j > and parity P, there_is a degenerate doublet 



of states with spin-parity J p = (j ± \) p . Following FalkE^a, I discuss the cases of 
integral and half integral j separately. Hadronic states with integral j have odd 
fermion number and correspond to baryons; states with half-integral j have even 
fermion number and correspond to mesons. 

First consider the heavy baryons. In this case, the "brown muck" is an object 
with spin-parity j p that can be represented by a totally symmetric, traceless tensor 
A** 1 "' 1 ** subject to the transversality condition A Ml = 0. States are said to have 
"natural" parity if P = (— l)- 7 , and "unnatural" parity otherwise. The composite 
heavy baryon can be represented by the tensor wave function 

c" -" ///, .1" -"' . (179) 

Under a connected Lorentz transformation A, this object transforms as a spinor- 
tensor field 

^w-w _^ a^ 1 . . . Ag D{A) i\) ux - v i , (180) 

where D(A) = exp^—jU^a^) is the usual spinor representation of A. A heavy 
quark spin rotation A, on the other hand, acts only on Uh', hence 

jjj^H _> D(A) V> Ml -^ . (181) 

Here A is restricted to spatial rotations (in the rest frame). The infinitesimal form of 
D(A) was considered in ( |117| ). The simplest but important case j p = + corresponds 
to the ground-state Ag-baryon with total spin-parity J p = | .It can be represented 
by a spinor u\. Since the light degrees of freedom are in a configuration of total 
spin zero, the spin of the baryon is carried by the heavy quark, and the spinor u\ 
coincides with the heavy quark spinor. Hence 

tpK = u A (v , s) = u h (v , s) . (182) 

For j > 0, the object ijj^ 1 —^ does not transform irreducibly under the Lorentz 
group, but is a linear combination of two components with total spin j ± \. These 
correspond to a degenerate doublet of physical states, which only differ in the orien- 
tation of the heavy quark spin relative to the angular momentum of the light degrees 
of freedom. The nonrelativistic quark model suggests that one should identify the 
states with j p = 1 + with the Sq (J p = | + ) and Sq (J p = | + ) baryons. In the 
quark model, these states contain a heavy quark and a light vector diquark with 
no orbital angular momentum. Unlike the Ag-baryons, they have unnatural parity. 

This i mpl ies that decays between Aq and Yiq* must be described by parity-odd form 
factorsHS 
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Next consider the heavy mesons. I shall only discuss the case j 1 



i- 

2 



m 



detail. Since quarks and antiquarks have opposite intrinsic parity, the corresponding 
physical states with "natural" parity are the ground-state pseudoscalar ( J p = 0~) 
and vector ( J p = 1~) mesons. As before, the heavy quark is represented by a spinor 
u h (v, s) subject to the condition $u h (v, s) = u h (v, s). The light degrees of freedom 
as a whole transform under the Lorentz group as an antiquark moving at velocity 
v. They are described by an antifermion spinor ve(v, s') satisfying vg(v, s') f = 
—vg(v,s'). The ground-state mesons can be represented by the composite object 
ip = Uh vg, which is a 4 x 4 Dirac matrix with two spinor indices, one for the heavy 
quark and one for the "brown muck". Under a connected Lorentz transformation 
A, the meson wave function ip transforms as 



i) -> D(A)iPD- 1 (A) , 
whereas under a heavy quark spin rotation A 

V> -> D(A) ip . 



(183) 



(184) 



The composite ip represents a linear combination of the physical pseudoscalar and 
vector meson states. It is easiest to identify these states in the rest frame, where Uh 
has only upper components, whereas ve. has only lower components. The nonvan- 
ishing components of ijj are thus contained in a 2 x 2 matrix, which can be written 
as a linear combination of the identity / and the Pauli matrices a\ Let me choose 
the quantization axis in 3-direction and work with the rest-frame spinor basis 
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Then a basis of states is: 



( U + ) = 
( U - ^ ) = 

V2(U) = 




x/2 V° 



(186) 



Let me furthermore define two transverse polarization vectors e± and a longitudinal 
polarization vector e 3 by 



4 = -^(0,l,±i,0), e^ = (0,0,0,l). 



(187) 
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It is th en obvious to identify the pseudoscalar (P) and vector (V) meson states 



a: 



204205 



1 + 7° 

P(V = 0) = --^ L 1B, 

V(v = 0,e)= ±±^. (188) 

The second state in (|186| ) has longitudinal polarization, whereas the last two states 
have transverse polarization. 

To get familiar with this representation, consider the action of the spin operator 
S on A matrix representation of the components S l in the rest frame is X* = 
§757°7 l j an d the action of the operator S* on the meson wave function is "S 1 ip = 

■?/>]• Using this, one finds 

f2p = S 3p = q ^ 

S 2 V(e) = 2V(e), 
S 3 V(e ± ) = ±V(e ± ), 

S 3 U(e 3 ) = 0, (189) 

which shows that P has total spin zero, and V has total spin one. Next consider the 
action of the heavy quark spin operator S. It has the same matrix representation 
as S, but only acts on the heavy quark spinor in ip: S l ip = S l ip, with S l = E\ It 
follows that 

S 3 P = \ V(e 3 ) , 
S 3 V(e 3 )= 1 -P 1 

S i V(e ± ) = ±±V(e ± ), (190) 

in accordance with the spin assignments for the heavy quark in ( |186|) . 

In a general frame, the tensor wave functions in ( 188 ) can be readily generalized 



in a Lorentz-covariant way by replacing 7 with f. The covariant representation of 
states can be used to determine in a very efficient way the structure of hadronic 
matrix elements in the effective theory. The goal is to find a minimal form factor 
decomposition consistent with Lorentz covariance, parity invariance of the strong 
interactions, and heavy quark symmetry. The flavour symmetry is manifest when 
one uses mass-independent wave functions. The correct transformation properties 
under the spin symmetry are guaranteed when one collects a spin-doublet of states 
into a single object. In case of the ground-state pseudoscalar and vector mesons, 
for instance, one introduces a covariant tensor wave function M.(v) that represents 
both P{v) and V(v, e) by 

ka/ \ _ 1+^ f — 75 ; pseudoscalar meson, (191) 
2 I ( ; vector meson. 
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It has the important property A4(v) = P + M.(v) P_, where P± = |(1 ± f). This is 
often used to simplify expressions. 

Consider now a transition between two heavy mesons, M(v) — > M'(v'), mediated 
by a renormalized effective current operator T ft^ , which changes a heavy quark Q 
into another heavy quark Q' . According to the Feynman rules of HQET, the "heavy 
quark part" of the decay amplitude is simply proportional to u' h T u^, interactions 
of the heavy quarks with gluons do not modify the Dirac structure of T. Since 
the heavy quark spinors are part of the tensor wave functions associated with the 
hadron states, it follows that the amplitude must be proportional to M. {y') Y JA(v). 
This is a Dirac matrix with two indices representing the light degrees of freedom. 
Since the total matrix element is a scalar, these indices must be contracted with 
those of a matrix E. Hence one may write 

(M'{v')\h' v ,Th v \M{v)) = Ti{E{v,v',fi)M\v')TM(v)}. (192) 

The matrix E contains all long-distance dynamics. It is a most complicated object, 
only constrained by the symmetries of the effective theory. Heavy quark symmetry 
requires that it be independent of the spins and masses of the heavy quarks, as well 
as of the Dirac structure of the current. Hence, E can only be a function of the 
meson velocities and of the renormalization scale /i. Lorentz covariance and parity 
invariance imply that H must transform as a scalar with even parity. This allows 
the decomposition 

E{v, v', fi) = Si + E 2 f + ~ 3 f + H 4 f f , (193) 

with coefficients Sj — Ei(v- v', fi). But using the projection properties of the tensor 
wave functions, one finds that under the trace 

E(v, v', fi) -> Hi - E 2 - ~ 3 + E 4 = -£(v ■ v', /x) . (194) 

Therefore!, 

(M'(v')\h' v ,Th v \M{v)) = -^vv',fx)Ti{M\v')TM(v)}. (195) 

The sign is chosen such that the universal form factor £(v ■ v', y) coincides with the 
Isgur-Wise function, which is the single form factor that describes semileptonic weak 
decay processes of heavy mesons in the infinite quark mass limit. Equation ( |195|) 
summarizes in a compact way the results derived in Section 3.4. Using the explicit 
form of the meson wave functions, one can readily recover the relations ( |128|) , ( |129| ). 
and (|138| ). The only new feature is that the Isgur-Wise function depends on the 



renormalization scale /i. This is necessary to compensate the scale dependence of 
the Wilson coefficients, which multiply the renormalized current operators in the 
short-distance expansion. 
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3.9. Meson Decay Form Factors 



One of the most important applications of heavy quark symmetry is to derive re- 
lations between the form factors parametrizing the exclusive weak decays B — > D Iv 
and B — ► D*lv. A detailed theoretical understanding of these processes is neces- 
sary for a reliable determination of the element | V c b\ of the Kobayashi-Maskawa 
matrix. Let me start by introducing a set of six hadronic form factors hi(w), which 
parametrize the relevant meson matrix elements of the flavour-changing vector and 
axial vector currents = cj^b and N 1 = 07^75 h: 

(D(v')\ V» \B{v)) = h + {w) (v + v'Y + h-(w) (v - v'f , 
{D*{v',e)\V^\B{v))=ih v {w)e^elv' a vp, (196) 
(D*(v', e)| A" \B(v)) = h Al (w) (w + 1) e*>*-[h M {w) v" + h A3 (w) v^]e*-v . 

Here w = v ■ v' is the velocity transfer of the mesons, and e denotes the polarization 
vector of the D*-meson. At leading order in the I/tuq expansion, one can derive 
expressions for these form factors by using the short-distance expansion (|169j ) of the 
flavour-changing currents, as well as the tensor formalism outlined above. According 
to ( |177p , the fj, dependence of the Wilson coefficients of any bilinear heavy quark 
current can be factorized into a universal function Ki^(w,fi), which is normalized 
at zero recoil. The \i dependence of this function has to cancel against that of 
the Isgur-Wise function. One can use this fact to define a renormalization-group 
invariant Isgur-Wise form factor as 

£ren(w) = £(u>, A*) ^hh(^, A 4 ) , Cren(l) = 1 • (197) 

Neglecting terms of order 1/rriQ, one then obtains^ 

(M'(v')\J^\M(v)) = -Un{w) Cf ] H Tr { M\v') Ti M(v)}, (198) 

i=i 

where the Dirac matrices have been defined in ( |170| ) and ( |171j ). It is now straight- 
forward to evaluate the traces to find 



h+(w) 
h-(w) 
h v (w) 



h A Aw) 



h M (w) 



h As {w) 



C x {w) + 
to + lr 



w + 1 



C 2 (w) + C 3 (w) 



C 2 (w) - C 3 (u?)J £ ren («0 



(199) 
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This is the correct generalization of ( |129| ) and ( |138[ ) in the presence of short-distance 
corrections. The fact that, to leading order in 1/rriQ, the meson form factors are 
givenin terms of a single universal function ^ ren (^) was the discovery of Isgur and 
Wise». Short-distance QCD corrections affect the form factors in a calculable way. 
Their effects are contained in the various combinations of short-distance coefficients, 
which can be evaluated using the numerical results given in Table [I]. 



3.10. Power Corrections and Luke's Theorem 



Using the covariant tensor formalism and the short-distance expansions of the 
effective Lagrangian and currents beyond the leading order in l/mq, one can inves- 
tigate in a systematic way the structure of power corrections to the relations derived 
in the previous section. I have given a simple (tree-level) example of the structure 
of power corrections in ( |143| ). For the more complicated case of meson weak decay 
form factors, the analysis at order l/mq was performed by Luke". Lately short- 
distance corrections were included to all orders in perturbation theory r 03 l r 08 l. Falk 
and myself have also analysed the structure of corrections for both meson 

and baryon weak decay form factor st53. I shall not discuss these rather technical 
issues in detail, but nevertheless give you the main results. 

Luke has shown that, for transitions between two heavy ground-state (pseu- 
doscalar or vector) mesons, the 1/toq corrections can be parametrized by a set of 
four additional universal functions of the velocity transfer w. The most important 
outcome of his analysis concerns the zero recoil limit, where an analogue of the 
Ademollo-Gatto theorem^ can be proved. This is Luke's theorem!!, which states 
that the matrix elements describing the leading 1/mg corrections to meson decay 
amplitudes vanish at zero recoil. As a consequence, in the limit v = v ' there are 
no terms of order I/tuq in the hadronic matrix elements in (|196| ). This theorem is 
independent of the struc 
in perturbation theory 



187 



HE 



20, 



oi the Wilson coefficients and thus valid to all orders 



20!) 



There is considerable confusion in the literature about the implications of this 
result. It is often claimed that the theorem would protect any meson decay rate, or 
even all form factors that are normalized in the spin- flavour symmetry limit, from 
first-order power corrections at zero recoil. However, these claims are erroneous. 
The reason is simple but somewhat subtle. Luke's theorem only protects the form 
factors h + and Iia 1 in ( p.96[ ), since all the others are multiplied by kinematic factors 
which vanish for v = v'. In fact, an explicit calculation shows that the 1/itiq 
corrections to hy, h,A 2 , and Ka 3 do not vanish at zero recoil. The fact that these 
functions are kinematically suppressed does not imply that they could not contribute 
to physical decay rates. This is often overlooked. Consider, as an example, the 
process B — > D i v in the limit of vanishing lepton mass. By angular momentum 
conservation, the two pseudoscalar mesons must be in a relative p-wave in order to 
match the helicities of the lepton pair. The amplitude is proportional to the velocity 
\vd\ of the D-meson in the 5-meson rest frame, which leads to a factor (w 2 — 1) in 
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the decay rate-Jji such a situation, form factors that are kinematically suppressed 



can contribute^. Indeed, the B — > D tv decay rate is proportional to 

(200) 



[w 2 - I) 3 / 2 



m B - m D 

h +{w) ■ h-{w) 

m B + m D 



Both form factors, h + and contribute with similar strength to the rate near 
w = 1, although h- is kinematically suppressed in (|196j ). Consequently, the decay 
rate at zero recoil does receive corrections of order 1/mq. The situation is different 
for B — > D*£ v transitions^ 2 !. Because the vector meson has spin one, the decay can 
proceed in an s-wave, and there is no helicity suppression near zero recoil. One finds 
that close to w — 1 the decay rate is proportional to (w 2 — l) 1 / 2 |/iAi(^)| 2 - Since 
the form factor h Al is protected by Luke's theorem, these transitions are ideally 
suited for a precision measurement of | V^ J fr om an extrapolation of the momentum 
spectrum of the D*-meson to zero recoilEj'tHa. There, the normalization of the decay 
rate is known in a model-independent way up to corrections of order 1/toq: 

h Al {l) = VA + 5 1/m 2 + ... , (201) 

where 77^4 ~ 0.99 contains the short-distance corrections [cf. (|178|) 1- 

One expects higher-order power corrections to be of order 5 x / m 2 ~ (Aqcd/^c) 2 ~ 
3%, but of course such a naive estimate could be too optimistic. For a precision 
measurement of | V c b \ , it is important to know the structure of such corrections in 
more detail. Although in principle straightforward, the analysis of I/tuq correc- 
tions in HQET is a tedious enterprisefilll. Three classes of corrections have to be 
distinguished: matrix elements of local dimension-five current operators, "mixed" 
corrections resulting from the combination of corrections to the current and to the 
Lagrangian, and corrections from one or two insertions of operators from the effec- 
tive Lagrangian into matrix elements of the leading-order currents. Within these 
classes, one can distinguish corrections proportional to l/wi 2 , ^/m 2 , or l/mhm c . 
More than thirty new universal functions are necessary to parametrize the second- 
order power corrections to meson form factors. When radiative corrections are 
neglected, eleven combinations of these functions contribute to the hadronic form 
factors hi(w). The general results greatly simplify at zero recoil, however. There 
the equation of motion and the Ward identities of the effective theory can be used 
to prove that matrix elements of local dimension-five current operators, as well as 
matrix elements of time-ordered products containing a dimension-four current and 
an insertion from the effective Lagrangian, can be expressed in terms of only two 
parameters Ai and A2, which are related to the I/tuq corrections to the physical 
meson masses. Moreover, the conservation of the flavour-diagonal vector current 
in the full theory forces certain combinations of the universal functions to vanish 
at zero recoil. The consequence is that whenever a form factor is protected by 
Luke's theorem, the structure of second-order power corrections at zero recoil be- 



comes rather simple. The quantity #i/ m 2 in (|201) is of this type. A careful estimate 
based on some mild model assumptions givesll^feHl —3% < 5i/ m a < —1%. One can 
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combine this result with (|201| ) to obtain one of the most important, and certainly 
most precise predictions of HQET: 

h Al (l) = 0.97 ±0.04. (202) 

3.11. Properties of the Isgur-Wise Function 

The establishment of heavy quark symmetry as an exact limit of the strong 
interactions enables one to derive approximate relations between decay amplitudes, 
and normalization conditions for certain form factors, which are similar to the 
relations and normalization conditions that can be derived for Goldstone-boson 
scattering amplitudes from the low-energy theorems of current algebra. HQET 
provides the framework for a systematic investigation of the corrections to the limit 
of an exact spin- flavour symmetry. The output of such a model-independent analysis 
is a short-distance expansion of decay amplitudes, in which the dependence on the 
heavy quark masses is explicit. At each order in the 1/mQ expansion, the long- 
distance physics is parametrized by a minimal set of universal form factors, which 
are independent of the heavy quark masses. As presented so far, the analysis is 
completely model-independent. Since hadronic decay processes are of a genuine 
nonperturbative nature, however, it is clear that predictions that can be made based 
on symmetries only are limited. In particular, very little can be said on general 
grounds about the properties of the form factors of the effective theory. But there 
is a lot of information contained in these functions. Much like the hadron structure 
functions probed in deep-inelastic scattering, they are fundamental nonperturbative 
quantities in QCD, which describe the properties of the light degrees of freedom in 
the background of the colour field provided by the heavy quarks. Since a static 
colour source is the most direct way to probe the strong interactions of quarks and 
gluons at large distances, a theoretical understanding of the universal form factors 
would not only enlarge the predictive power of HQET, but would also teach us in 
a very direct way about the nonperturbative nature of the strong interactions. 

The leading-order Isgur-Wise function £,(w) plays a central role in the descrip- 
tion of the weak decays of heavy mesons. It contains the long-distance physics 
associated with the strong interactions of the light degrees of freedom and cannot 
be calculated from first principles. Nevertheless, some important properties of this 
function can be derived on general grounds, such as its normalization at zero recoil, 
which is a consequence of current conservation. According to ( |128|) , the Isgur-Wise 
function is the elastic form factor of a ground-state heavy meson in the limit where 
power corrections are negligible. As such, £(u>) must be a monotonically decreasing 
function of the velocity transfer w = v ■ v ', which is analytic in the cut w-plane 
with a branch point at w = — 1, corresponding to the threshold q 2 = 4:mQ for heavy 
quark pair production. However, being obtained from a limiting procedure, the 
Isgur-Wise function can have stronger singularities than the physical elastic form 
factor. In fact, the short-distance corrections contained in the function K^(w,fi) 
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lead to an essential singularity at w = — 1 in the renormalized Isgur-Wise function 
denned in (|T9^ ). 



When using a phenomenological parametrization of the universal form factor, 
one should incorporate the above properties. Some legitimate forms suggested in 
the literature are: 

2 r , o , w - 1 



2 r 
£bswM = — — r exp ^ - (2g - 1) 

U> + 1 I W + 1 

£isgw(™) = exp { - £ 2 (w - 1 



2 W 



epoieH=(^ TT ) • (203) 



The first function is the form factor derived by Rieckert and mysel£i3 from an 
analysis of the BSW mo deli, the second one corresponds to the ISGW mo 
and the third one is a pole-type ansatz. Of particular interest is the behaviour 
of the Isgur-Wise function close to zero recoil, which is determined by the slope 
parameter g 2 > defined by £'(1) = — g 2 , so that 

£{w) = 1- g 2 (w-l) + 0[(w-l) 2 ]. (204) 

It is important to realize that the kinematic region accessible in semileptonic decays 
is small (1 < w < 1.6). As long as g 2 is the same, different functional forms of £,(w) 
will give similar results. A precise knowledge of the slope parameter would thus 
basically determine the Isgur-Wise function in the physical region. 

Bjorken has shown that g 2 is related to the form factors of transitions of a 
ground-state heavy meson into excited states, in which the light degrees of freedom 
carry quantum numbers j p = | + or | + , by a sum rule which is an expression of 
quark-hadron duality: In the infinite mass limit, the inclusive sum of the proba- 
bilities for decays into hadronic states is equal to the probability for the free quark 
transiti on If one normalizes the latter probability to unity, the sum rule has the 
formS 



212 



213 



i = ^{ieHI 2 + Eie (/) HI 2 } 

+ {w-i) {2E I^HI 2 + (v + i) 2 E k$MI : 

+ 0[(w-l) 2 ], (205) 



where /, m, n label the radial excitations of states with the same spin-parity quantum 
numbers. The sums are understood in a generalized sense as sums over discrete 
states and integrals over continuum states. The terms on the right-hand side of the 
sum rule in the first line correspond to transitions into states with "brown muck" 
quantum numbers j p = \ . The ground state gives a contribution proportional to 
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the Isgur-Wise function, and excited states contribute proportionally to analogous 
functions £®(w). Because at zero recoil these states must be orthogonal to the 
ground state, it follows that £^(1) = 0, and one can conclude that the corresponding 
contributions to ( |205| ) are of order (w — l) 2 . The contributions in the second line 



correspond to transitions into states with j p = | or |\ Because of the change in 
parity, these are p-wave transitions. The amplitudes are proportional to the velocity 
\vf \ = (w 2 — l) 1 / 2 of the final state in the rest frame of the initial state, which explains 
the suppression factor (w — 1) in the decay probabilities. The functions Tj(w) are the 
analogues of the Isgur-Wise function for these transitions. Transitions into excited 
states with quantum numbers other than the above proceed via higher partial waves 
and are suppressed by at least a factor (w — 1) 2 . 

For w = 1, eq. ( |205| ) reduces to the normalization condition for the Isgur-Wise 
function. The Bjorken sum rule is obtained by expanding in powers of (w — 1) 
and keeping terms of first order. T aki ng into account the definition of the slope 



parameter g 2 in (|2U4|) , one finds thatcMELi 



Q 2 = \ + H Irf^HI 2 + 2£ \r^(w)\ 2 > \ . (206) 



Notice that the lower bound is due to the prefactor \(w + 1) of the first term in ( |205| ) 
and is of purely kinematic origin. In the analogous sum rule for Aq-baryons, this 
factor is absent, and consequently the slope p aramet er of the baryon Isgur-Wise 



function is only subject to the trivial constraintl 210 lr 14 l g 2 > 0. 

Based on various model calculations, there is a general belief that the con- 
tributions of excited states in the Bjorken sum rule are sizeable, and that g 2 is 
substantially larger than 1/4. For instance, Blok and Shifman have estimated the 
contributions of the lowest-lying excited states to (|206| ) using QCD sum rules and 
find that§l! 0.35 < g 2 < 1.15. The experimental observation that semileptonic B- 
decays into excited D**-mesons have a large branching ratio of abouti3 2.5% gives 
further support to the importance of such contributions. 

Voloshin has derived another sum rule involving the form factors for transitions 
into excited states, which is the analogue of the "optical sum rule" for the dipole 



scattering of light in atomic physics. It readsEIZ 



mM mQ = J>M |r$Ml a + 2J>$ \4i(w)\ 2 , (207) 

^ m n 

where Ej are the excitation energies relative to the mass rriM of the ground-state 
heavy meson. The important point is that one can combine this relation with the 
Bjorken sum rule to obtain an upper bound for the slope parameter g 2 : 

(208) 

where E m i n denotes the minimum excitation energy. In the quark model, one expects^ 
that E m i n ~ mjf — mg, and one may use this as an estimate to obtain g 2 < 0.75. 

§ Strictly speaking, the lowest excited "state" contributing to the sum rule is D + %, which has an 



67 



The above discussion of the sum rules ignores renormalization effects. However, 
both the slope parameter g 2 in ( |206| ) and the heavy quark mass rriQ in (|207| ) are 
renormalization-scheme dependent quantities. Although there exist some q\ 
tive ideas of how to account for the \i dependence of the various parameters 
there is currently no known way to include renormalization effects quantitatively. 
One should therefore consider the bounds on g 2 as somewhat uncertain. To account 
for this, I relax the upper bound derived from the Voloshin sum rule and conclude 
that 



0.25 < g 2 < 1.0 



(209) 



where it is expected that the actual value is close to the upper bound. Recently, 
de Rafael and Taron claimecLip have derived an upper bound g 2 < 0.48 from 
general analyticity arguments^. If true, this had severely constrained the form 
of the Isgur-Wise function near zero recoil. It took quite some time to become 
clear what went wrong with their derivation: The effects of resonances b_elpw the 



threshold for heavy meson pair production invalidate the argument^ ~Ki£fl. It is 
possible to derive a new bound, which takes into account the known properties of 
the T-statesBH3. However, it is too loose to be of any phenomenological relevance, 
and one is thus left with the sum rule result (p09|) . 
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Figure 16: QCD sum rule prediction for the Isgur-Wise 
function in the kinematic region accessible in semileptonic 
decays^. The dashed lines indicate the bounds on the slope 
at v • v' = 1 derived by Bjorken and Voloshin [see (|209|) 1. 

Any sensible calculation of the Isgur-Wise function has to respect these bounds. 



As an example, I show in Fig. 16 the result of a QCD sum rule analysis of the 



excitation energy spectrum with a threshold at m T . However, one expects that this spectrum is 
broad, so that this contribution will not invalidate the upper bound for g 2 derived here. 
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renormalized universal form factorEH E3. Similar predictions have been obtained for 
the functions that describe the 1/mg corrections to the meson form factors^. 



3.12. Model- Independent Determination of \ V c b\ 

One of the most important results of HQET is the prediction ( |202|) of the nor- 
malization of the hadronic form factor tiA 1 at zero recoil. It can be used to obtain 
a model-independent measurement of the element | V c b\ of the Kobayashi-Maskaj^a 
matrix. The semileptonic decay B — > D*£i? is ideally suited for this purpose^. 
Experimentally, this is a particularly clean mode, since the reconstruction of the 
.D*-meson mass provides a powerful rejection against background. From the the- 
oretical point of view, it is ideal since the decay rate at zero recoil is protected 
by Luke's theorem against first-order power corrections in 1/mq. In terms of the 
hadronic form factors denned in ( |196| ), one finds 

1 dT(B^D*£u) G 2 F \V cb \\ s2 o l7 ,. t2 

hm 7 == -L— 1 = -Jd_^L {mB _ m ^f mh , {hAi . (210 ) 

At zero recoil, the normalization of Ha x is known in a model-independent way with 
an accuracy of 4%. Ideally, then, one can extract | V c b\ with a theoretical uncertainty 
of well below 10% from an extrapolation of the spectrum to w = 1. 

Presently, the proposal to measure | V c b\ close to zero recoil poses quite a chal- 
lenge to the experimentalists. First, there is the fact that the decay rate vanishes at 
zero recoil because of phase space. Therefore the statistics gets worse as one tries 
to measure close to w = 1. However, I do not believe that this will be an important 
limitation of the method. The phase-space suppression is proportional to \/w 2 — 1 
and is in fact a rather mild one. When going from the endpoint wj max — 1.5 down 
to w = 1.05, the change in the statistical error in | V c b\ 2 due to the variation of the 
phase-space factor is not even a factor of 2. A more serious problem is related to 
the fact that, for experiments working on the T(4s) resonance, the zero- recoil limit 
corresponds to a situation where both the B- and the D*-mesons are approximately 
at rest in the laboratory. Then the pion in the subsequent decay D* — ► D n is 
very soft and can hardly be detected. Thus, the present experiments have to make 
cuts which disfavour the zero recoil region, leading to large systematic uncertainties 
for values of w smaller than about 1.15. This second problem would be absent at 
an asymmetric 5-factory, where the rest frame of the parent 5-meson is boosted 
relative to the laboratory frame. 

In view of these difficulties, one presently has to rely on an extrapolation over a 
wide range in w to obtain a measurement of | V^|. In general, the differential decay 
rate can be written asH'llH 

— - — — (m B - m D *) m 6 D , rj A Vw 2 - 1 (w + 1) 



dw 487T 3 

4w 1 — 2wr + r 2 l . T r l9 - 9 
1 + ^^^ Zr- Wcb\ 2 Z 2 (w), (211) 



w + 1 (1 — r) 
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where tja = 0.99 is the short-distance correction to the form factor }ia 1 {w) at zero 
recoil, and r = mf)./mg. Equation ( |211| ) is written in such a way that the deviations 



from the heavy quark symmetry limit are absorbed into the form factor £(w), which 
in the absence of symmetry-breaking corrections would be the Isgur-Wise function. 
Since everything except | V c b\ and is known, a measurement of the differential 
decay rate is equivalent to a measurement of the product | V c b\ £(w). However, theory 
predicts the normalization of at zero recoil: 

= Va 1 M 1 ) = 1 + S 1/m 2 = 0.98 ± 0.04 , (212) 

where the uncertainty comes from power corrections of order l/m%. Using this 

information, | V c b\ and £(w) can be obtained separately from a measurement of the 
differential decay rate. 

I have applied this strategy for the first time-iH to the combined sample!^ of the 
data on B° -> D* + £ v decays collected until 1989 by the ARGUS and CLEO col- 
laborations. For the extrapolation to zero recoil, I used the parametrizations given 
in (|203| ) for the function £(w), treating its slope at zero recoil as a free parameter. 
The result obtained for | V c b\ was 

T B \V2 



iKfcl — = 0.039 ±0.006. (213) 
VI. 5 ps/ 

Since this original analysis the data have changed. In particular, the branching ratio 
for D* + — > D° 7r + has increased fromS 55% toUH 68%. This lowers the decay rate 
for B° — ► D* + £u, and correspondingly decreases | by 10%. However, the new 
data recently reported by the ARGUSlll and CLEOE§ collaborations give a larger 
branching ratio than the old data, indicating that further changes in the analysis 
must have taken place. It is thus not possible to simply rescale the result for | V c b\ 
given above. 



In Fig. [T7|, I show the new ARGUS dataEi^ for the product | V c b\ £,(w). From an 



unconstrained fit using again the parametrizations in ( |203| ) , the following value is 
obtained: 



iKil — — = 0.049 ±0.008. (214) 
VI. 5 ps/ 

However, the fit gives very large values for the slope parameter g 2 , between 1.9 and 
2.3. Since the slope of agrees with the slope of the Isgur-Wise function up 
to power corrections of order Aqcd/^c, I believe that such large values cannot be 



tolerated in view of the Voloshin sum rule. In Fig. |17|, I therefore show a fit to the 
data which uses the pole ansatz in ( |203j ) with the constraint that g 2 < 1. This leads 
to a significantly smaller value: 



VJ — — = 0.037 ±0.006. (215) 
VI. 5 ps/ 
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Figure 17: ARGUS dataN for the product 



function of the recoil w 



v - v , assuming r^o 



Vcb\ £(w) as a 
: 1.5 ps. Using 



the normalization condition ( |212p , \V c b\ is obtained from an 
extrapolation to w = 1. The fit curve is explained in the text. 



Very recently, the CLEO collaboration has reported new results for the recoil 
spectrum with higher statistics^. They have applied tight cuts in order to reduce 
the systematic errors. These data are shown in Fig. [18]. The curves correspond to 
the results of a fit using the various functions in ( |203| ) . Without any constraint on 
the slope at zero recoil, one obtains 

I V cb \ (^^-\ 12 = 0.039 ±0.006, (216) 
VI. 5 ps/ 

together with values 1.0 ± 0.4 < g 2 < 1.2 ± 0.7. It is reassuring that the result for 
the slope parameter is in agreement with the bound derived from the Voloshin sum 
rule. Because of the careful analysis of systematic uncertainties, I consider these 
results to be the best numbers for | V c b\ and g 2 that are currently available. 

The model-independent determination of | V c b\ is one of the many applications 
of heavy quark symmetry that have been explored over the last years. I hope that 
the above presentation gives you a flavour of the potential of the HQET formalism. 
For a more comprehensive discussion, I refer to my review papered. 

3.13. Exercises 



• Go through the derivation of the effective Lagrangian of HQET and prove 
eqs. (JTOl-Q. 
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Figure 18: CLEO data@ for the product \ V cb \£(w) as a 
function of the recoil w. 



Derive the Feynman rules depicted in Fig. 12 from the effective Lagrangian 

O- 

Convince yourself that ( |138| ) is the correct covariant expression for the pseudo- 
scalar-to- vector transition matrix element. 

Derive the one-loop expression given in ( |149| ) for the wave-function renormal- 
ization constant Z^. 

By performing the appropriate Dirac traces in (|198|) , derive the expressions 
given in ( |199| ) for the meson form factors hi(w). 

Derive the Bjorken sum rule (|2U6p by taking the limit w — > 1 in (I 
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Concluding Remarks 

In these lectures, I have presented an introduction to the theory and phenomenol- 
ogy of heavy quark masses and mixing angles, as well as to the exciting new field 
of heavy quark spin-flavour symmetry. The synthesis of these topics is interesting 
in itself, in that it combines rather mature subjects with one of the most active 
fields of research in particle physics. As far as quark masses and mixing angles are 
concerned, the theoretical concepts have been set many years ago. Yet new theo- 
retical tools have to be developed to determine these parameters more and more 
accurately. Heavy quark effective field theory, on the other hand, has received broad 
attention in recent years only, although some ideas of a spin-flavour symmetry for 
heavy quarks had been around for much longer. These new developments start to 
have an impact on the precision with which one can determine some elements of 
the Kobayashi-Maskawa matrix. I have discussed this in detail for the extraction 
of | V c b\ from the semileptonic decays B — > D*iv. 

Besides reviewing the status of the theoretical developments, my purpose was to 
convince you that -B-physics is a rich and diverse area of research. Presently, this 
field is characterized by a very fruitful interplay between theory and experiments, 
which has led to many significant discoveries and developments on both sides. Heavy 
quark physics has the potential to determine important parameters of the flavour 
sector of the electroweak theory, and at the same time provides stringent tests of 
the standard model. On the other hand, weak decays of hadrons containing a heavy 
quark are an ideal laboratory to study the nature of nonperturbative phenomena 
in QCD. 

When I presented these lectures in the summer of 1993, the future of -B-physics 
was uncertain. Today, the prospects for further significant developments in this 
field look rather promising. With the approval of the first asymmetric _B-factory at 
SLAC, and with ongoing .B-physics programs at the existing facilities at Cornell, 
Fermilab, and LEP, there are clearly .Beautiful times ahead of us! 
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